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ABSTRACT 


This report presents a simple and rapid method for predicting the 
nonlinear response and stiffener separation of postbuckled, flat, 
composite, shear panels. The method uses modular, “strength-of- 
materials-type” models. The disbonding failure is hypothesized as being 
due to a stress concentration in the surface layer of the skin immediately 
under the toe of the attached flange. The highly local character of this 
stress, which renders a finite element analysis impractical, enables a 
simple analysis to be developed. The maximum principal tensile stress in 
the skin surface layer under the toe is determined, and failure is said to 
occur when this stress reaches the mean transverse tensile strength of 
the layer. 

The analysis consists of a number of closed-form equations that 
have been programmed into a preliminary design code called SNAPPS 
[Speedy Nonlinear Analysis of Postbuckled Panels in Shear], which rapidly 
predicts the panel’s postbuckling response and failure load. SNAPPS was 
applied to three test panels with widely different geometries, laminates 
and stiffnesses. At the test failure loads, the predictions of maximum 
principal tensile stress at the toe vary from 8% below to 20% above the 
transverse tensile strength. The predicted failure loads range from 18% 
below to 25% above the test failure loads. 

Using the B-basis allowable stress, which is 62% of the mean 
transverse tensile strength, the predicted allowable load for each panel 
was consistently found to be about half the test failure load. In view of 
the scatter experienced in matrix-dominated modes of failure in 
composite materials, these are reasonable, but not excessively 
conservative predictions. 

A new method was developed for determining experimental buckling 
loads for shear panels. The method uses the "diagonal tension factor", k, 
for which a closed-form expression was derived. An unambiguous 
estimate of the initial buckling load can be obtained by plotting k against 
load, and extrapolating back to the load at which k is zero. 

The report also gives design recommendations for reducing the 
maximum principal tensile stress at the toe, thereby increasing the 
panel’s load carrying capacity. Results from an approximate but accurate 
parametric analysis imply this is best accomplished by decreasing the 
skin thickness, increasing the flange width, and by decreasing the bending 
stiffness of the combined stack of the skin and attached flange. 
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SECTION 1 


INTRODUCTION 

There is a need for a simple, rapid prediction of stiffener separation 
in postbuckled, flat, composite shear panels in which the stiffeners are 
bonded to the skin. Analysis and test evidence point to the failure being 
caused by a highly localized stress in the skin immediately under the toe 
of the attached flange of the stiffener. To locate and quantify this stress, 
a finite element model would have to be exquisitely detailed and may not 
be practical in real terms, particularly if it has to be re-done at each 
design change. Fortunately, the extreme localization of the stress enables 
the development of a simplified method of analysis. 

1.1 Methodology 

The analysis described in this report began with the premise that 
the physical behavior of a highly buckled, stiffened shear panel could be 
described by simple equations culled from standard texts. The challenge 
was to get a physical understanding of the panel behavior as it progressed 
into the postbuckled regime. Two key aids to this understanding came 
with a series of panels tested by Grumman [Visconti, 1988], followed by a 
detailed finite element analysis [Sobel and Sharp, 1 994] of one test panel, 
denoted as Cl. That analysis used the STAGS code [Almroth, et al., 1982] 
to perform a buckling analysis and then a nonlinear analysis up to the test 
failure load of fifteen times initial buckling. 

The analysis developed here consists of a number of closed-form 
equations that can easily be used in a "hand analysis". For expediency, 
they have been programmed using the True Basic language into a code 
called SNAPPS [Speedy Nonlinear Analysis of Eostbuckled Panels in Shear], 
which rapidly predicts the postbuckling response and failure load of the 
panel. The model used is similar to that of Tsai [1983], where strips of 
the skin and stiffener parallel to the diagonal tension field are idealized 
as beams and tie-rods. We depart from Tsai’s analysis in the derivation 
of the out-of-plane displacement, in the idealization of the combined skin 
and flange, and in the failure criterion. For the latter, we compute a 
highly localized stress in the skin immediately under the toe of the flange 
and relate it to the transverse tensile strength of the skin ply, while Tsai 
uses an overall pull-off load which is geometry-dependent and has to be 
determined empirically for each configuration. 
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• Development Phases 

The development of SNAPPS can be divided into three phases: 

The first phase drew heavily on the methods of Kuhn [1952], Wagner 
(described in Kuhn [1952]), and Ranalli and Bunce [1977], for the analysis 
of stiffened aluminum panels in which the stiffeners were riveted to the 
panel. Ranalli and Bunce's analysis focused on overall pull-off loads 
between the skin and flange. This concept was carried over to a test 
program for stiffened composite panels, test coupons, and elements, as 
described by Visconti [1988]. Their analysis, when modified to include the 
twisting moment between the skin and flange in addition to the pull-off 
load [Sharp and Sobel, 1989] gave satisfactory agreement with test 
failure loads, but required element tests specific to the test panel, and 
these showed adhesive failures. In contrast, the panel tests revealed 
intra-ply failures of the skin immediately under the toe of the flange, as 
evidenced by the pull-out of fibers from the surface layer of the skin. 
This observation led to the realization that the flatwise tensile strength 
in the laminate was the key to the failure. 

The second phase began with an analysis of the through-the- 
thickness tensile stress in the skin immediately under the toe (Sharp 
[1989] and Sobel [1990]). The model used was developed during Grumman’s 
CTSA program [Cacho-Negrete, 1978], and consists of two beams, one 
representing the skin and the other the attached flange, joined by an 
elastic foundation. The analysis revealed peel stress concentrations 
under the heel and toe of the flange. These stresses die away within short 
distances, and do not significantly interact with each other, enabling a 
great simplification in the equations. In this second phase, the stress 
state in the postbuckled skin was assumed to be given by Kuhn's empirical 
theory of incomplete diagonal tension, and the width of the buckle was 
determined empirically from moire' fringe patterns. 

The third phase took account of the immense amount of data 
available from the STAGS analysis of panel Cl, Sobel [1990], together 
with a more detailed interpretation of strain gauge readings and moire' 
fringe patterns for the three test panels. From these, it was found that a 
reasonably accurate estimate of the stress resultants in the flat panel 
after buckling could be obtained by assuming (1) a diagonal tension angle 
of 45° and, (2) that the compressive stress resultant reaches a value equal 
to the buckling shear flow and does not increase further. The buckle width 
could now be determined analytically by isolating one buckle and imposing 
equilibrium on the out-of-plane components of the membrane stress 


1-2 



resultants. All of the empiricism in the earlier analysis was thereby 
removed. 

• Failure Criterion 

It is hypothesized that the failure is located in the skin layers 
immediately under the toe of the flange, and is caused by the maximum 
principal tensile stress in the layers (Sharp [1989] and Sobel [1990]). 
This stress is a combination of three stresses, all of which peak at the 
toe, in the outermost skin layer on the stiffener side of the panel whose 
fibers are normal to the diagonal tension fold. They are (1) an extremely 
localized peel stress, (2) an interlaminar shear stress, and (3) an inplane 
tensile stress in the direction of the diagonal tension fold. The 
combination of these three stresses makes the toe location critical: at the 
heel, the peel stress is higher, but the combination of the three stresses 
gives a lower principal stress there. The failure criterion adopted here 
equates the maximum principal tensile stress based on these three 
components to the transverse tensile strength of the skin layers. 

• Assessment of Method 

An assessment of the accuracy of the simplified method is achieved 
by applying the SNAPPS code, described in Appendix M, to the three panels 
manufactured and tested by Grumman [Visconti, 1988]. The predictions 
are compared to those of a nonlinear STAGS finite element analysis [Sobel 
& Sharp, 1994], and to strain gauge and moire" fringe readings taken during 
the tests. SNAPPS predictions for the out-of-plane displacement, 
wavelengths, and skin stress resultants were found to agree well with 
both test measurements and STAGS results. While the intensely localized 
state of stress at the toe cannot be directly compared to the STAGS 
results, nor can it be determined from the test results, relating its 
predicted maximum principal tensile value to the transverse tensile 
strength of the skin layers provides estimations of the panel strength that 
are in reasonably good agreement with the test failure loads. 

1.2 Structure of Report 

This report consists of a number of sections and appendices, with 
details of the derivations of the governing equations used in SNAPPS being 
relegated to the appendices. The governing equations and computational 
procedure are summarized in Section 2, and coded in appendix M. Section 3 
compares SNAPPS predictions for the three test panels with test 
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measurements and STAGS results, and Section 4 gives conclusions, 
lessons learned and recommendations. 

The key elements of SNAPPS involve calculation of the following: 

• The stress resultant state in the postbuckled skin (Appendix C), 

• The buckle width and length (Appendix D), 

• The maximum out-of-plane displacement (Appendix E), 

• The “tie-rod” analysis of a strip of skin (Appendix F), 

• The skin bending moment at the toe of the flange from 
compatibility of rotation between the tie-rod model and the 
combined skin and flange model at the toe of the flange 
(Appendix G), 

• The maximum shear stress at the toe (Appendix H), 

• The tensile peel stress between the skin and the toe of the flange, 
based on a simplified "coupled beams on an elastic foundation" 
analysis (Appendix I), 

• The inplane transverse tensile stress in the critical ply, which is 
the outermost skin ply whose fibers are normal to the diagonal 
tension fold (Appendix K), 

• The combining of the peel stress, the maximum shear stress and 
the transverse tensile stress under the toe to obtain the maximum 
principal tensile stress in the surface layer of the skin 
(Appendix L), 

• Failure, which is identified as value of the applied shear flow at 
which the maximum principal tensile stress reaches the 
transverse tensile strength of the critical ply (Appendix L). 
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SECTION 2 


DESCRIPTION OF ANALYSIS & COMPUTATIONAL PROCEDURE 
2.1 Description of Analysis 

The analysis described here is based on four major insights gained 
from the test results and the STAGS analysis: 

(1) That the angle of the diagonal tension fold remained almost 
constant at between 45° and 48° to the normal to the stiffener. 

(2) That the stress resultants in the skin after buckling could be 
simply described in the terms of the applied shear flow and the 
buckling shear flow without recourse to empiricism. 

(3) That the buckles grew steadily in the skin as the load increased 
above buckling by narrowing the width of the buckle. This 
observation identified the aspect ratio and maximum out-of-plane 
displacement of the buckle as key parameters in the analysis. 

(4) That the failure was localized under the toe of the attached 
flange of the stiffener, enabling a simplification of the CTSA 
[Cacho-Negrete 1978] analysis to be used to calculate the tensile 
peel stress in the skin under the toe. 

Based on these insights, a simplified computational procedure was 
developed that uses strength-of-materials-type models. The first step in 
the procedure is to calculate the buckle width and length as a function of 
the applied shear flow. Next, the maximum out-of-plane displacement is 
computed. From this, the bending moment and transverse shear force in the 
skin immediately under the toe of the stiffener flange are determined. With 
these internal loads at the toe region, the analysis proceeds to the 
calculation of three stresses, t xz , cr z and o x , all of which peak at the 
surface of the skin under the toe of the flange. They are defined in the xz 
plane, where x is in the direction of the diagonal tension fold and z is 
normal to the original plane of the skin. The shear stress, t xz , arises from 
the transfer of part of the diagonal tension load from the skin into the 
flange, c z is a tensile peel stress acting normal to the skin between the 
skin and the flange, and a x is a tensile stress in the outermost ply on the 
stiffener side of the panel whose fibers are normal to the direction of the 
diagonal tension fold. o x is induced by the membrane load and bending 
moment in the skin immediately under the toe. Combining these stresses 
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enables calculation of the maximum principal tensile stress in the critical 
ply. Comparing this to the transverse failure stress of the ply allows a 
prediction of failure to be made. 

2.2 Governing Equations & Computational Procedure 

In what follows, only the main equations are given. Detailed 
derivations of the equations are in the referenced appendices, and a flow 
chart of the various steps in the analysis procedure is given in Figure 1. 

Note that the input to SNAPPS requires the membrane and bending 
stiffnesses of the skin, flange and the combined skin and flange (where they 
are bonded together to form a single stack). These have to be obtained from 
the layer properties and stacking sequences, using an independent computer 
code. We used MACLAMINATE [Flanagan, 1991], but any similar capability 
will serve. 

• Skin Stress Resultants 

Next, the principal stress resultants in the buckled skin, except under 
the toe of the flange, are determined from Appendix C as follows: 

Ni = 2q-q cr Eq. (7) of Appendix C 

N 2 = qcr Eq. (5) of Appendix C 

where Ni and N 2 are principal tensile stress resultants parallel and normal, 
respectively to the direction of the diagonal tension folds. Ni is positive 
when it is in tension and N 2 is positive when it is in compression. In these 
equations, q is the applied shear flow, and q cr is the shear flow at which 
the skin buckles (q and q C r are always positive in this report). The buckling 
shear flow q C r may be estimated from standard procedures for orthotropic 
plates, such as the Grumman Advanced Composites Structures Manual 
(which is based on Seydel [1933]) or Housner and Stein [1975]. If available, 
the test buckling load or a STAGS (or similar) prediction can be used. 

• Buckle Kinematics 

The length, L, of the diagonal tension fold, measured along the fold 
between the toes of the attached flanges of the stiffeners, see Figure 1 of 
Appendix D, is given by 


L = btt/cosa 


2 -2 



where 

btt - the toe-to-toe distance between the attached flanges of the 
stiffeners, measured normal to the stiffener direction, 
and 

a = the angle between the diagonal tension fold and the normal to the 
stiffener direction. It is assumed that a is constant at 45°. 


The width of the buckle, c in Figure 1 of Appendix D, is calculated 


from 



Eq. (7) of Appendix D 


The projected length of the buckle along the stiffener axis, Lb, is then 

Lk - — 2 — Eq. (8) of Appendix D 

u cos a 


• Maximum Out-of-Plane Displacement 

The amplitude of the out-of-plane displacement, 5 m ax> is next 
obtained from 

5max = Eq. (6) of Appendix E 

where Gt is the in-plane shear stiffness of the skin, referred to axes 
parallel and normal to the stiffener direction. 

• Loads in the Skin at the Toe of the Attached Flange 

The analysis now obtains the bending moment, Mt 0 e. and the transverse 
shear, Wtoe. in the skin immediately under the toe of the flange. To do this, 
a "tie-rod" model of the skin is connected to a beam model representing the 
combined flange and skin between the stiffener centerline and the toe of 
the flange. This is done by invoking rotational compatibility at the common 
ends of the models. Both models are of unit width and lie along the diagonal 
tension fold. 

The "tie-rod" model of the skin is described in Appendix F. The beam 
model of the combined skin and flange, which contains the total stacking 
sequence of the flange, a layer of adhesive and the skin, is described in 
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Appendix G. It is assumed to be cantilevered at the stiffener centerline, 
and loaded at it's tip by Wtoe and a moment given by (Mtoe - Nitoez) The 
second term, NitoeZ. represents the moment due to the membrane load in the 
unbuckled skin at the toe being offset from the neutral axis of the combined 
section by an amount, z, (see Figure 1 of appendix G) and is calculated from 


and 


where 


Nitoe = q 


Etfrts + t a +-^- + Etsr-f- 

2 / 2 . is. 

Etfr + Et sr 2 


Eq. (10) of Appendix C 


Eq. (6) of Appendix G 


t s = thickness of the skin 
tf = thickness of the flange 

t a = thickness of the adhesive layer between the skin and flange 
Etsr * membrane stiffness of the skin 
Etfr = membrane stiffness of the flange. 


Et sr and Etf r are measured in the direction of the diagonal tension fold (as 
denoted by the subscript r). 

Mtoe is given by 


ft 5max | -j 

M»e = j-!= L. 

[ A, - te/L . 
I Ni 

and Wtoe by 

Wtoe = ®-(Ni Smax 

where 



and 

Lfr = Lf/cosa, where Lf ■ the width of the attached flange of the 
stiffener, measured normal to the stiffener direction 

D S r = the bending stiffness of the skin about the inplane normal to the 
diagonal tension fold 


LfrNt 

2 Dtotr 


Lfr q z 
Dtotr 


Lf, 


JL 


Lfr 


Dtotr 2 L Dtotr^ 


Eq. (10) of Appendix G 


+ Mtoe) 


Eq. (15) of Appendix G 
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Dtotr = the bending stiffness of the combined skin and flange (where 
they are bonded together to form a single stack). D to tr is 
measured about the inplane normal to diagonal tension fold. 

• Calculation of Stresses in the Skin at the Toe 

The analysis now proceeds to the calculation of the three stresses 
which combine to give the maximum principal tensile stress at the surface 
of the skin under the toe. The stresses are defined in the xz plane, where x 
is in the direction of the diagonal tension fold and z is normal to the 
original plane of the skin, and are as follows: 

The shear stress x X z at the interface between the skin and the flange 
is calculated first. It arises from the transfer of part of the membrane 
load from the skin into the flange and peaks in the surface ply under the toe. 
From Appendix H, the maximum shear stress, which occurs under the toe, i s 
given by 


'txz — Nitoe K— 

Etfr + Et sr 


where 


K = -a / <|)xz jrr — ~ — 
V \Etfr Etsr 


and 


<j>xz = 


1/ t 




2\Gi3/fiange 'G'adhesive 



Eq. (5) of Appendix H 
Eq. (4) of Appendix H 

Eq. (2) of Appendix J 


In the equation for 4> xz , G denotes the shear modulus for the isotropic 
adhesive and core materials. For the tape layers, it is simplest to use G 13 
for all layers regardless of orientation, while for fabric layers G 13 = Gi 2 - 

The peel stress o z is a tensile stress acting normal to the skin face 
between the skin and the flange. This stress has peaks at the toe and heel 
of the flange, but we are concerned only with the peak at the toe in this 
analysis. As described in Appendix I, o 2 can be conveniently subdivided into 
two parts: the first, a z i, is due directly to the moment Mtoe and the 
transverse shear Wtoe applied from the skin; and the second, a Z 2 , is caused 
by the shear along the interface between the skin and flange acting at half 
the laminate thickness away from the midplane of each laminate. The two 
parts of the peel stress are given by 
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Gz- * + O z 2 


where 


^*(TT^) (Wt - + PMtoe) 

and 


Eq. (21) of Appendix I 


0 Z 2 = Cn tf fe(Fi + F2) 


where 


R _ / foz Dsr + Dfr 
\ 4 D sr Dfr 


il / 4 


Q n _ Nitoe K Etfr 

2 4>z Etfr + Etsr 


Fi = 


1 - mlPfr K _. ■ 1 


i 2 2 p' 

iZP 

Dfr + L»sr 


♦j 


p_ = — P / P sr/Dfr ~ ts/tf \ 
K \ Dfr + Dsr / 

0z 


m = 


D S r Df r 


^ Dsr K 4 + foz |2 ^ s ^f) 

K 4 + 4 p 4 j 


Eq. (22) of Appendix I 


and 


<t>z = 


i|J-) + 

2 \ E 3 /flange 


1 


i) + (j_) + lit.) 

E 'adhesive ’E 3 /face 2'E' C ore 


Eq. (1) of Appendix J 


In these equations, Df r is the bending stiffness of the flange, 
measured about the normal to diagonal tension fold, t denotes laminate 
thickness, E is Young’s modulus for the isotropic adhesive and core 
materials, and E3 is the out-of-plane modulus for the flange and skin layers 
(by assuming transverse isotropy for tape layers, we may take E3 = E2, 
where E2 is the layer transverse modulus 


The tensile stress o x in the skin face is induced by the membrane 
stress resultant and bending moment in the skin immediately under the toe 
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of the flange. It acts in the direction of the diagonal tension fold, in the 
"critical ply". The critical ply is defined by the following criteria: it has 

to be on the stiffener side of the skin; it has to be put into tension, along 
the diagonal tension fold, by the transverse shear Wtoe and the moment Mtoe 
and; it's fibers are normal to the diagonal tension fold, i.e. g x acts normal 
to the fibers, and in the plane of, the critical ply. To be conservative, by 
using the highest possible value of c x , the critical ply is assumed to be at 
the surface of the skin. The force and bending moment are analyzed i n 
Appendix G. The derivation of g x is given in Appendix K, and results in the 
following equation for a x : 




E 2 

(l "V 1 2 V 2 1) 


(l-V 12 )(l +V S r) 


q 

Etsr 


+ t§^ Mtoe 
2 Dsr 


Eq. (7) of Appendix K 


where 

v S r = the Poisson's ratio of the skin in the xy axis system (v S r = v xy 
is defined according to v xy E y = Vy X E x , where E x and E y are the 
engineering constants for the skin laminate), 


Ei, E 2 , and Vi 2 are the in-plane Young's moduli and Poisson’s ratio 
for the face ply of the skin. (The subscript 1, when affixed to a 
material property, pertains to the fiber direction, and v i 2 is defined 
according to v 2 i Ei = v -|2 E 2 ). 


• Maximum Principal Tensile Stress 

The three stress components x xz , o z and o x are combined to give the 
following expression for the maximum principal tensile stress in the 
surface ply immediately under the toe: 

a mpt = + s\j + (^z) 2 Eq. (1 ) of Appendix L 


• Failure Criterion 

Examination of photomicrographs of ply cross sections shows that for 
graphite/epoxy the placing of the fibers within such a section is random, 
and it is impossible to discern the orientation of the photographed section 
within the ply without other, external, clues. We may therefore assume 
transverse isotropy, and that the tensile strength at anv orientation in the 

xz plane is equal to the transverse tensile strength F$> u of the layer. 
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Failure is assumed to occur when the maximum principal tensile 
stress in the skin immediately under the toe reaches the transverse tensile 
strength of the layer, i.e., when 

G mpt = F 2 U Eq. (2) of Appendix L 

Failure of this layer constitutes failure of the joint, because cracking of 
the matrix allows fibers to be pulled out of the skin surface. This type of 
failure was observed on all of the panels tested. 

The following values used in the analysis of the test panels are based 
on statistical analysis of 19 coupon test results for transverse tension of 
IM6/3501 -6 graphite/epoxy tape, in the room temperature, ambient, 
moisture (“dry”) condition [Shyprykevich, 1988]: 

F* 2 U - Mean strength = 7150 psi 
Standard Deviation = 1180 psi 
"B-basis" allowable stress = 4460 psi 

where the B-basis allowable strength is such that at least 90% of the 
transverse tensile strengths are expected to exceed the B-basis allowable 
value, with a confidence of 95%. 
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Specify New Value of q 


SNAPPS INPUT : (Appendix M) 

Specify Geometry, Material Properties, Stiffness, q cr 


Specify the Applied Shear Flow, q, and Compute the Following: 


Stress Resultants N 1( N 2 (Appendix C) 


Wavewidth, c and Wavelength, (Appendix D) 


Maximum Out-Of-Plane Displacement, S max (Appendix E) 


Skin Bending Moment at Toe of Flange, M toe and Transverse Shear Force 
at Toe, W toe (Appendix G) 


Foundation Moduli, <j> 2 , <}) xz (Appendix J) 


Stress Components At Toe 

Interlaminar Shear Stress 

^xz 

(Appendix H) 

Peel Stress 
(Appendix 1) 

Tensile Stress in Critical 
Ply, in DT Direction, O x 

(Appendix K) 


Maximum Principal Tension Stress at toe, CJ mpt (Appendix L) 
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SECTION 3 


RESULTS AND ASSESSMENT OF METHODOLOGY 

The computational procedure described in Section 2 has been 
programmed, using the True Basic language, into a code called SNAPPS, 
(Appendix M), which has been applied to three, flat, stiffened, composite 
panels (A1, B1 and Cl) tested at Grumman, as described in Appendix A. 
The current section presents results pertaining to the SNAPPS 
methodology, and assesses the accuracy of the SNAPPS predictions by 
comparing them to the panel test results and STAGS predictions for panel 
Cl [Sobel and Sharp, 1994], 

3.1 Assessment of Basic Premises 

The SNAPPS analysis is built around two basic premises: that the 
angle of diagonal tension, a, is 45° (for flat panels), and that the 
compressive stress resultant normal to the diagonal tension fold reaches 
a value equal to the buckling shear flow, qcr> but does not increase beyond 
that value as the panel is loaded into the postbuckling regime. The first 
assumption implies that the principal stress tensile and compressive 
resultants, Ni , and N 2 , line up with the diagonal tension folds at all load 
levels, with Ni acting in the direction of the fold and N 2 normal to it, and 
the second assumption stems from the use of Wagner’s model of the 
buckled panel, see Kuhn [1952], Inplane equilibrium then dictates the 
value of Ni, which is derived in Appendix C. The validity of these 
assumptions is assessed by comparing SNAPPS values for a, Ni and N 2 
with test and STAGS results for the three panels. This is done in Figures 
1 through 4, with the caption of Figure 1 defining the location of the 
strain gauges, and how the gauge measurements at these locations were 
averaged to obtain the test values of oc, Ni and N 2 . (The values of the test 
failure loads and q C r appearing in the abscissas in the figures are given in 
Appendix A). 

Figure 1 compares the a = 45° assumption with measured and STAGS 
results. The measured and predicted values of a all lie between 45° and 
48° degrees from the normal to the stiffener, thereby justifying this 
assumption. 

Figures 2, 3 and 4 plot the principal stress resultants for the three 
panels against q/q C r- Test results are given together with predictions 
based on the SNAPPS model and STAGS (for Cl only). Also included are 
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results from Kuhn’s [1952] model based on an empirical expression for the 
diagonal tension factor, k, obtained from tests on aluminum panels. The 
figures reveal that the simplified SNAPPS approach gives stress 

resultants that are in reasonable agreement with the measured and STAGS 
results. 

From Figures 1 through 4 it is concluded that the basic premises 
used in the SNAPPS analysis correlate well with the test and STAGS 
results. SNAPPS results stemming from these basic premises are 
compared next with test and STAGS for both global (kinematic) and local 
(stress) response variables. 

3.2 Assessment of Global Predictions (Kinematic Variables) 

SNAPPS computes the buckle (half) wavelength projected along the 
stiffener direction, Lb, (Appendix D) and the maximum out-of-plane 
displacement, 5 m ax (Appendix E). Normalized values of these kinematic 
variables are compared next with test and STAGS results to provide a 
"global" check on the SNAPPS methodology. 

Figure 5 gives the variation of Lb with the dimensionless load level 
q/q cr . Lb, is normalized with respect to the toe-to-toe distance, btt, 
measured normal to the stiffener direction. This distance is shown in 
Figure 1 of Appendix D, and the SNAPPS expression for Lb/ btt is given by 
Eq. (10) of that appendix. Included in the figure are wavelengths measured 
from moire' photographs taken during the loading to failure for panels A1 
and Cl, and wavelengths predicted by the STAGS analysis for Cl, and 
Timoshenko’s [1961] “end point (q/qcr =1)” buckling solution through 
which the SNAPPS prediction was forced to pass (see Appendix D). The 
figure demonstrates that the simple SNAPPS model used to determine the 
wavelength gives results that agree quite well with the measured and 
STAGS results. 

Figure 6 shows the variation with q/q cr of 5 ma x» which is normalized 
with respect to LVycr, where L is the toe-to-toe distance between 
flanges, measured along the length of the diagonal tension fold (see Figure 
1 of Appendix E), and Ycr is the critical shear strain corresponding to q C r- 
The SNAPPS expression for LYy C r is given by Eq. (8) of Appendix E. The 
test results in the figure are based on moire' fringe patterns photographed 
at various load levels during the tests of panels A1 and Cl (panel B1 was 
the earliest tested and was not instrumented for moire' fringe 
measurements). The test set-up for the Shadow Moire’ Method of obtaining 
out-of-plane displacements in the buckled panel is described in Visconti 
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[1988]. The SNAPPS predictions based on the simple model of Appendix E 
somewhat overpredict 5 m ax- but this discrepancy is judged acceptable 
within the framework of a simplified, preliminary design, analysis tool. 

Taken together, Figures 1 through 6 engender a sense of confidence 
that the basic premises underlying the SNAPPS procedure are physically 
realistic, and are valid even at postbuckling load levels of fifteen times 
initial buckling (which corresponds to the failure load level for panel Cl). 

3.3 Assessment of Local Predictions (S tress Variables) 

The final check involves using SNAPPS to predict the maximum 
principal tensile stress, o mp t, immediately under the toe of the flange, 
and comparing this stress to the transverse tensile failure stress of the 
skin ply. This stress is highly localized and cannot be verified, either by 
the tests or STAGS results. However, if we get reasonable predictions for 
the failure loads for the three panels, then we may reasonably conclude 
that the methodology is valid, at least within the range of parameters for 
the tested panels, for the purpose of providing a preliminary design tool. 
The comparisons are made in two ways, the results of which are 
summarized in Tables 1 and 2. 

Corresponding to the test failure load for each of the three tested 
panels, Table 1 compares o mp t with the mean transverse tensile strength, 

F 2 U , which is taken as the failure criterion, as described Section 2 and 
Appendix L. For the IM6/3501-6 Graphite/Epoxy material in the room 

temperature, ambient moisture content condition, F^ = 7150 psi (see 
Appendix L). Table 1 shows that <J mpt differs from the mean strength by an 
average of 12.9%. It can be shown that the calculated values of <T mpt are 

approximately within one standard deviation from the mean value of F^. 

The SNAPPS code can also be used to predict the failure shear flow. 
To do this, the code is run with a succession of values of the applied shear 
flow until the predicted maximum principal tensile stress at the toe 
reaches the mean strength of 7150 psi. The results so obtained are 
displayed in Table 2, from which it is seen that the predicted failures 
differ from the test failures by an average of 18.3% (based on absolute 
differences, using algebraic differences gives -2%). This is a higher 
percent difference than the 12.9% obtained in the previous table, because 
o m pt varies nonlinearly with q. 
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Table 2 also gives SNAPPS predictions of the allowable load for each 
panel based on the "B-basis" allowable of 4460 psi, (which is 62% of the 
mean failure stress of 7150 psi, see Appendix L). It is observed that the 
predicted allowable loads for the three panels are consistently about half 
the test failure load. These predicted allowable loads are reasonably but 
not excessively conservative, in light of the nonlinear nature of the 
problem, and the variability in the matrix-dominated strength. 

Tables 1 and 2 give SNAPPS results for a mpt for two load levels, one 
corresponding to the test failure load for each panel, and the other 
corresponding to the load level at which SNAPPS predicts failure to occur. 
The SNAPPS predictions for o mpt for a range of load levels, are displayed 
in Figures 7 to 9 for the three panels, which also serve to summarize the 
foregoing results of Tables 1 and 2. The SNAPPS predictions for the 
maximum principal tensile stress under the toe are seen to vary 
nonlinearly with load, thereby demonstrating that the simple models used 
in SNAPPS capture the basic geometric nonlinear behavior of the 
postbuckled shear panels. 



TABLE 1 SNAPPS PREDICTIONS FOR THE MAXIMUM PRINCIPAL 
TENSILE STRESS AT THE TEST FAILURE LOAD FOR EACH OF THE 
THREE TESTED PANELS. 


Panel 

Predicted Max Prin 



Stress at toe, psi 


A 1 

8597 

20.2 

B 1 

7893 

10.4 

Cl 

6580 

-8.0 


Avg % Diff of Absolute 

12.9 


Values of o m pt Relative to 



Mean Strength 



TABLE 2 SNAPPS FAILURE AND ALLOWABLE LOAD PREDICTIONS 
FOR EACH OF THE THREE TESTED PANELS. 


Panel 

Test Failure 

Predicted 

% Diff of 

Allowable 

Ratio of 


Load, ppi 

Failure 
Load, ppi 

Predicted 
Load w.r.t 
Test 

Shear 
Flow (B- 
basis), 

Allowable 
Load to Tesl 
Failure Load 




Fail ure 
Load 

ppi 


■M 

1403 

1149 

-18.1 

791 

.56 

■3P 

838 

735 

-12.3 

451 

.54 

SEE* 

962 

1199 

24.6 

468 

.49 


Avg % Diff of 
Absolute 
Values of 
Predicted 
Failure Load 
w.r.t to Test 
Failure Load 

18.3 
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Note: The Angle of Diagonal Tension Is Computed From Strain Gage Values at the 

Gage Locations In Each Bay. These Values are Then Averaged at the Following Locations: 

• For Cl (see Inset): Center of the Left and Right Bays. 

• For A & B1: 1/4 and Center Points In Center Bay, and 1/4 Point In End Bay. 
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FIGURE 2 VARIATION OF PRINCIPAL STRESS RESULTANTS WITH 
APPLIED SHEAR LOAD FOR PANEL Cl (q =63.5 ppi FROM STAGS) 
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APPLIED SHEAR LOAD FOR PANEL A1 (q =634 ppi FROM TEST) 
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FIGURE 5 VARIATION OF AXIAL WAVELENGTH OF BUCKLE 
WITH APPLIED SHEAR FLOW IN SKIN 
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FIGURE 6 VARIATION OF MAXIMUM OUT-OF-PLANE 
DISPLACEMENT OF BUCKLE WITH APPLIED SHEAR FLOW IN SKIN 
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& PREDICTED FAILURE LOAD. PANEL A1. 

Note the nonlinear variation of predicted stress with load. 
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FIGURE 8 MAXIMUM PRINCIPAL TENSILE STRESS AT TOE, 
& PREDICTED FAILURE LOAD. PANEL B1. 
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FIGURE 9 MAXIMUM PRINCIPAL TENSILE STRESS AT TOE, 
& PREDICTED FAILURE LOAD. PANEL Cl. 








SECTION 4 


CLOSING REMARKS 


4.1 Summary and Conclusions 

We have attempted, using a simple "strength-of-materials" 
approach, to develop a preliminary design tool for rapidly analyzing the 
very complex nonlinear problem of the postbuckling behavior and failure of 
flat, stiffened, composite shear panels. 

The analysis is based on the premise that the physical behavior of a 
highly buckled, stiffened shear panel can be described by simple, closed- 
form, equations using the insights gained from a series of panels tested 
by Grumman, and the subsequent correlation of the test results with 
STAGS nonlinear finite element results for one of the panels. 

The analysis method has been programmed, using the "True Basic" 
language, into a code called SNAPPS ("Speedy Nonlinear Analysis of 
Eostbuckled Eanels in Shear"), which is a simple and fast preliminary 
design tool for predicting failure and postbuckling behavior of composite 
shear panels with bonded stiffeners. The code is written in a modular 
form so that it may be easily modified by the user if new test data or 
physical insight makes modification desirable. Each module of the code i s 
the subject of an Appendix in this report. 

Confirmation of the validity of the basic premises of the method has 
been gained by comparing its predictions with test results for the three 
panels, A1, B1, and Cl, and with the nonlinear finite element analysis 
results obtained from STAGS for panel Cl only. We are encouraged that 
the SNAPPS predictions for the skin’s diagonal tension stress resultants 
and out-of-plane deflections (i.e., the large-scale behavior of the panel) 
agree well with test results for these panels of widely different 
geometries and stiffnesses, and with the STAGS predictions. 

However, the failure criterion adopted - that failure occurs when 
the maximum principal tensile stress in the skin under the toe of the 
stiffener flange exceeds the transverse tensile strength of the surface 
plies - cannot be confirmed, either by test measurements or by finite 
element analysis. Indirectly, it is judged by comparing the test failure 
loads of the A1, B1 and Cl panels with SNAPPS predictions based on the 


4 -1 



mean strength of 19 transverse tensile test coupons. The predictions of 
failure for the two-stiffener panels, A1 and B1, are 18% and 12% below 
the test failures respectively, i.e. on the conservative side. For panel Cl , 
the predicted failure load is about 25% above the test failure load, i.e., 
unconservative. The mean (algebraic) error is 2%. We contend that the 
failure criterion is physically realistic, and that the failure mode is 
correctly identified, and feel this conclusion is confirmed by Minquet and 
O’Brien [1996] in their analysis of stiffener pull-off specimens. It is 
worthy of note that their analysis of the maximum principal tensile 
stress at failure for four different specimens gave values ranging from 
6600 psi to 7300 psi for the same IM6/3501-6 graphite/epoxy material 
used in our tested panels and for which we used an average transverse 
tensile strength of 7150 psi. Never the less, we must admit that our data 
base is too small for comfort because of the large amount of scatter, 
although this is typical of matrix dominated failures. 

The B-basis allowable transverse tensile strength is 4460 psi, 
much lower than the mean value of 7150 psi (allowable = 62% of mean), 
reflecting the wide scatter of a matrix-dominated failure in 

graphite/epoxy, which is made worse by the small sample size. The 

prediction of allowable shear flows for all three panels is consistent, at 
56%, 54% and 49% of the test failure loads for A1, B1 and Cl respectively. 
We feel that these are reasonably, but not excessively, conservative in 
view of the matrix-dominated nature of the failure, and the small data 
base. 


SNAPPS requires much more work for its complete validation. Its 
predictions need to be compared with test results for other panels. A 
larger, statistically significant, data base of transverse tensile strengths 
needs to be acquired, and applied to a wider range of test panels. 
Meanwhile, SNAPPS should be useful in performing trade studies, with the 
understanding that any design obtained from the trade studies should be 
verified by test and/or a detailed analysis. 

4.2 Limitations 

The method was developed for flat , stiffened, composite panels 
loaded in shear only. The skin is assumed to be a balanced and symmetric 
laminate. 

There must be a sufficient number of stiffeners so that the aspect 
ratio of a skin-bay between the toes of the stiffener flanges is large 
enough to allow a number of buckles to develop, ensuring that end effects 
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do not dominate the buckle pattern, and that the out-of-plane deformation 
of the diagonal tension fold is symmetric about the midlength of the fold, 
as is assumed in the SNAPPS analysis (the STAGS results for panel Cl 
show that such symmetry is not the case for this single stiffener panel — 
see Figure 2 of Appendix G). 

The stiffeners are assumed to have sufficient bending stiffness to 
prevent the diagonal tension folds from progressing across the stiffener 
centerlines. The user is alerted that the analysis assumes zero rotation 
of the skin at the stiffener centerline in the diagonal tension direction 
(see Figure 1 of Appendix G). This is correct for stiffeners made up of 
back-to-back channels or angles, and may reasonably apply to sturdy 
"hat"-sections. The applicability of the code to unsymmetric sections 
such as single angles or Z-section stiffeners has not been verified, nor 
has it been applied to non-sandwich composite skins, or metallic panels 
with bonded-on stiffeners. 

The Wagner model [see Kuhn 1952] for the principal stress 
resultants in the diagonal tension folds has been used in this analysis 
because it gave good agreement with the measured and STAGS skin stress 
resultants, and it enabled the removal of undesirable empiricism from the 
analysis. This model should be checked against other test results. 

Further, because of the assumption that the principal stress 
resultants line up with the diagonal tension folds, and that the angle of 
diagonal tension is 45°, the method cannot be applied to curved panels. 

It is to be recalled that the allowable and mean transverse tensile 
strength values for the graphite/epoxy face ply material for the tested 
panel’s are based on a limited data base (see Appendix L). Thus, the values 
of the SNAPPS predictions for the allowable and failure loads may suffer 
accordingly. 

4.3 Lessons Learned & Design Implications 

Particular care needs to be taken with sign conventions, stacking 
sequence "clocks" on drawings, and load directions with respect to the 
structural axes (parallel and normal to stiffeners). Changes in stacking 
sequence or load direction relative to the structural axes on the drawing, 
resulting from manufacturing errors or changes in load spectrum, may 
have significant effects on initial buckling and failure. For thin skins, the 
initial buckling load is markedly effected by the sign of the shear load, 


4 -3 



and by orthotropy and/or anisotropy in the skin, as is well illustrated in 
Figure 17 of Nemeth [1997], 

While the stiffener must be sturdy enough to prevent the 

diagonal tension folds from progressing across the stiffener centerlines, 

the resistance to stiffener disbonding can be improved by reducing the 

moment in the skin immediately under the toe of the flange, thereby 

lowering the value of the maximum principal tensile stress. The moment 

is lower when the entire attached flange is flexible, but, because the 
stress components which make up the maximum principal tensile stress 
die away rapidly from the toe, the same effect may be achieved by 

tapering the flange locally at the toe, instead of the square edge we 

tested. Minquet and O'Brien [1996] have shown that trimming the flange 

edge to a 20° wedge angle at the toe gives a significant improvement in 
the load-carrying capability of the skin-stiffener joint. The taper may 
also be achieved by dropping internal plies in the stacking sequence of the 
flange, but this may cause problems in tooling because a terminated ply 
has to be prevented from moving during the cure cycle and because a resin 
pocket is created at it's end. 


4.4 Recommendations 

To further verify SNAPPS and to suggest possible refinements to it, 
SNAPPS should be applied to other test panels, particularly those without 
sandwich skins, those with many stiffeners and those tested in different 
fixtures. 

The models used in the various steps in the computational procedure 
are quite simple. Further development should concentrate on refinement 
of some of these models, particularly those associated with the transition 
region at the nose of the buckle where the diagonal tension field stress 
resultants in the buckled panel have to be transformed into a complex set 
of internal loads, both inplane and out-of-plane, at the toe of the flange. 
The refinement would be based on insights gained from a detailed finite 
element modeling of the area around the nose of the buckle where it meets 
the toe of the flange, and the model would have to include the shear 
deformation of the skin. The model would be bounded by adjacent nodal 
lines in the skin and should also include a portion of the attached flange 
bounded by the nodal lines, and the toe and heel of the stiffener. From the 
results of the model it would be possible to determine the redistribution 
of load and moment in the flange and the skin at the nose of the buckle, 
and thereby refine the current models used in SNAPPS. 



Although the above model would have to be super-detailed, it is still 
unlikely to accurately determine the concentrated stresses in the vicinity 
of the toe. In order to quantify the critical through-the-thickness 
maximum principal tensile stress in a one inch wide stiffener pull-off 
specimen, Minquet and O'Brien [1995 and 1996] used 2273 elements with 
13908 degrees of freedom. The size of the element was of the order of 
one-third of a ply thickness (i.e. about 0.002 inches) in the critical region. 
This was essentially a two-dimensional plane strain analysis where the 
location of the critical stress was known. Similar models of equal 
complexity were used by Wang, et al. [1994] and Li, et al. [1996]. Imagine, 
then, the complexity of the same finely detailed analysis in the present 
problem, where the analysis has to be three-dimensional. The critical 
stress is in the skin under the toe of the stiffener, but it's location could 
be anywhere along the length of the stiffener. Further, because of the 
changing wave-pattern in the postbuckled state, the location of the 
maximum stress will move with increasing load until, at some point and 
some load level, it reaches either the allowable stress for a design 
analysis, or the average strength for a failure analysis. Therefore, we 
should explore the possibility of using a code such as SUBLAM [Flanagan, 
1993] to investigate the stresses at the interface between the skin and 
the flange, using loads from the finite element model as boundary 
conditions in a global-local approach. 

The “diagonal tension factor method” (Appendix B) of obtaining the 
initial buckling shear load from experimental results needs additional 
verification. 

For panel Cl , which failed deep in the postbuckling range at a load 
about 15 times initial buckling, the predicted moment in the skin at toe, 
M toe , reached a maximum before the test failure load was reached and then 
decreased with increasing load (see Figure 3 of Appendix G). Further study 
is required to ascertain whether or not this behavior reflects physical 
reality or a limitation of the basic SNAPPS methodology. 

Finally, it is emphasized that only one failure mode is considered 
here, based on the observation that failure originated in the plies of the 
face sheet immediately under the toe of the flange for the three tested 
panels. Other failure locations and modes need to be investigated, such as 
those associated with the sandwich core for which the deformation of the 
skin at the nose of the buckle can be extreme at high postbuckling loads. 
Such an investigation is beyond the current scope of SNAPPS. However, 
the modular approach used in SNAPPS enables the user to estimate the 



transverse shear force and bending moment in the skin at the toe of the 
flange, and add more failure criteria to the code. 
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APPENDIX A 


DESCRIPTION OF TEST PANELS A1, B1, and Cl 

Comparisons are made throughout this report between SNAPPS 
predictions and test results for three shear panels tested by Grumman. 
This appendix briefly describes these panels: complete descriptions are 
given in the test report [Visconti, 1988], and a STAGS analysis report 
[Sobel and Sharp, 1994]. 

The panels (Figures 1 to 3) were tested in the room temperature 
"dry" condition (i.e. stabilized at the ambient moisture condition of the 
test laboratory). The test frame was designed to allow the diagonal 
tension field to be reacted by the compressive load in the stiffener(s), 
rather than in the test frame members. This is accomplished through the 
use of a "breather joint", as described below. The skin of each panel 
consists of two identical thin composite faces separated by SynCore to 
form a sandwich (Figure 2). SynCore is an epoxy resin filled with glass 
micro-balloons to reduce its density by approximately one half. The face 
material is IM6/3501-6 graphite/epoxy tape. The panels all have the 
same overall dimensions, so as to fit into the test frame, but the face 
layups and thicknesses differ, as does the thickness of the SynCore. 
Panels A1 and B1 have two stiffeners at approximately 11 inches spacing, 
but differ in the skin faces, the SynCore thickness and in the stiffener 
thickness. Panel Cl had a single stiffener, dividing the skin into two bays 
of approximately 17 inches width, and also differs from both A1 and B1 in 
its SynCore and stiffener thicknesses. The stiffeners (Figure 2) are built 
up from back-to-back channel sections of AS4/3501-6 graphite/epoxy 
fabric, with straps of the same layup as the channels. Because the 
longitudinal and transverse modulii are almost equal in a fabric, see Table 
1, the stiffener web and flanges are, as far as is possible with a wrapped 
shape, balanced and symmetric laminates. Warping in the curing process 
is thereby minimized. The stiffeners have the same overall dimensions, 
being formed on common tools, but are of different layups and thicknesses 
for the three panels. The laminate nominal thicknesses and layups are 
given in Tables 2 to 4. The stiffeners are cured first and then the "green" 
skin is cured and bonded to the stiffeners in one operation ("cocured") on a 
flat tool. A single layer of FM300 adhesive is used between the skin and 
the flange of the stiffener. 

The test frame is 38-in. -square fixture built up from heavy steel 
bars pinned at the frame corners. The test load P is applied through two 
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of these pins on opposite corners of the test frame, as shown in Figure 1 . 
Parallel to the stiffener(s), the thickened edges of the panel skin are 
sandwiched between the steel bars, and normal to the stiffener(s) the skin 
is attached to the steel bars through a breather joint. The breather joint 
(Figure 3) is intended to ensure that the diagonal tension field load in the 
buckled skin is reacted by compression loads in the stiffener(s) rather 
than in the test frame members. The breather joint consists of thin steel 
angles bolted to the thickened laminate at the edge of the skin, and these, 
in turn, are bolted to similar angles attached to the test frame. The joint 
acts as a bellows, capable of transmitting the applied shear load in the 
plane of the frame but not permitting significant axial load to develop 
across the joint. While the breather joint did induce axial compression in 
the stiffeners, it also contributed to the nonuniformity of shear stress 
near the panel edges, as illustrated in Figure 11 of Sobel and Sharp [1994]. 
Fortunately, the shear stress concentrations in the skin are well away 
from the disbond failures experienced between the skin and the attached 
flanges of the stiffener(s), and did not influence the test results. 

The panels were instrumented with strain gauges, using rosettes on 
both skin faces, and axial gauges on the stiffener web and flanges Layouts 
of the gauge positions for each panel are given in the test report. The 
moire fringe technique, described in Volume III of Visconti [1988], was 
used to obtain both wavelengths and out-of-plane displacements in the 
buckled skin for Panels A1 and Cl. Examination of a series of photographs 
of the moir£ fringes taken during the tests enabled the progression of the 
buckle waves to be recorded as the load increased. A major finding of the 
test program was the slow growth of new buckles progressing in from the 
ends of the panel, i.e. the buckle wavelength reduces continuously as the 
shear load is increased without sudden "snap-throughs". 

The mode of failure observed for all the test panels was not a 
failure of the bondline between the flange of the stiffener and the skin, 
but rather an intra-ply failure of the epoxy matrix of the skin immediately 
under the toe of the flange. This was evidenced by the pull-out of fibers 
from the surface layer of the skin. It was this observation that led to the 
realization that the flatwise tensile strength of the skin face material 
was the key to the failure. The test failure loads (from Visconti [1988]) 
and buckling loads (from Appendix B) are given in Table 5. 
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FIGURE 1 PLAN VIEW OF THE STIFFENER SIDE OF THE PANAL 
(ILLUSTRATED FOR THE SINGLE STIFFENER PANAL Cl ). 
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FIGURE 2 SKIN AND STIFFENER FOR PANELS A1 , B1 , AND Cl 
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(1) [Flanagan, 1988] 

(2) [Cox and Flanagan, 1997] 

(3) [Dastin, 1996] 

(4) E3=E2 assuming transverse isotropy for tape 

(5) V ]2 is defined according to V 12 E 2 =V 2 iEi 














TABLE 2 NOMINAL LAMINATE THICKNESSES (inches) 
OF TEST PANELS 


LAMINATE 

PANEL 

A 1 

B 1 

Cl 

EACH FACE 

.0336 

.0224 

.0224 

SYNCORE 

.050 

.040 

.020 

TOTAL SKIN 

.1172 

.0848 

.0648 

ADHESIVE 

.008 

.008 

.008 

STRAP 

.036 

.0288 

.0216 

WRAP 

.036 

.0288 

.0216 

TOTAL FLANGE (WRAP+STRAP) 


.0576 

.0432 

TOTAL WEB (WRAP+WRAP) 


.0576 

.0432 


TABLE 3 LAMINATE LAYUPS FOR TEST PANELS 


(a) Skin Lavup rM6/350l-6 Tape) 


LAMINATE 

PANEL i 


A 1 

B 1 

Cl 

Free Surface 




SKIN FACE 

45 




135 




0 

45 

45 


90 

135 

135 


135 

90 

90 


45 

0 

0 

SvnCore Surface 





Fabric) 






LAMINATE 

PANEL I 


A 1 

B 1 

Cl 

Stiffener Tool Face 




WRAP 

45 




90 

45 



0 

90 

45 


90 

0 

90 


135 

135 

___.U5.__. 

STRAP (inner & outer) 

135 

135 

45 


90 

0 

90 


0 

90 

135 


90 

45 



45 



MM 





NOTES: 1. 0° plies are parallel to stiffener direction (See Figure 1). 

2. The cross section and tool surfaces are shown in Figure 2. 
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TABLE 4 STACKING SEQUENCE FOR SKIN & ATTACHED 
FLANGE OF TEST PANELS. Stacking Sequences Are From 
Skin Tool Surface (See Figure 2), And Are Used In 
The Maclaminate Code To Obtain Stiffnesses For 
Input To Snapps. 


LAMINATE 

PANEL 

(Material) 

A 1 

B 1 

Cl 

FACE 

45 



(IM6/3501-6 Tape) 

135 




90 

45 

45 


0 

135 

135 


135 

90 

90 


45 

0 

0 


5 Plies 

4 Plies 

2 Plies 





FACE 

45 

0 

0 

(IM6/3501 -6 Tape) 

135 

90 

90 


0 

135 

135 


90 

45 

45 


135 




45 



ADHESIVE :l : lii r ; v ' 

m 


jpum 

(FM300) ... 

111 



STRAP 

45 



(AS4/3501 -6 Fabric) 

90 

45 



0 

90 

45 


90 

0 

90 


135 

135 

135 

WRAP 

'135 

135 

45 

(AS4/3501 -6 Fabric) 

90 

0 

90 


0 

90 

135 


90 

45 



45 




NOTE : 0° plies are parallel to stiffener direction (See Figure 1). 


TABLE 5 BUCKLING & TEST FAILURE LOADS FOR TEST PANELS 


LOAD 

PANEL 1 

A 1 

B 1 

Cl 

Buckling Load P C r, Kips 

37.5 

17.6 

4.13 

Failure Load Pfaiitest, Kips 

85.0 

50.0 

64.0 

Buckling Shear Flow q cr , ppi 

634. 

302. 

63.5 

Failure Shear Flow qtaiitest, PP> 

1403. 

838. 

962. 
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APPENDIX B 


BUCKLING SHEAR FLOW q cr FOR THREE TEST PANELS, AND 
DEVELOPMENT OF A NEW METHOD OF OBTAINING EXPERIMENTAL 
BUCKLING LOADS FOR FLAT SHEAR PANELS 

The SNAPPS code input requires the value of the buckling shear flow, 
q cr . This appendix describes methods used to obtain the buckling shear 
flow for the three panels A1, B1, and Cl (see Appendix A) analyzed here. 

B.1 Buckling Load for Panel C 1 

Sobel and Sharp [1994] performed a STAGS [Almroth, et al., 1982] 
analysis for one of the panels, panel Cl. The detailed finite element 
model includes the skin, stiffener, breather joint and loading frame. The 
STAGS results for Cl are as follows: 

• £LL P cr.STAGS = 41 3 Kips, q cr ,STAGS = 63.5 ppi (1 ) 

where P is the load applied at two opposite corners of the test frame (see 
Figure 1 of Appendix A). A new method, described next, was developed 
here for obtaining the buckling loads for the other two panels. 

B.2 Description of a Method of Obtaining Experimental Buckling Loads 
For Flat Shear Panels — Based on the Diagonal Tension Factor, k 

As is well known, the determination of buckling loads for flat panels 
from test results is notoriously difficult because the buckles initially 
develop slowly with little out-of-plane distortion and it is hard to 

recognize the load at which they first appear. This is particularly so i f 
the stiffeners are on one side of the skin only, so that bending out of the 
original flat plane of the skin occurs immediately upon application of 
load, and is aggravated if the as manufactured skin is not perfectly flat. 

Usually, the buckling load is estimated from plots of measured skin 

strains by finding the load at which the strains bifurcate or begin to grow. 
Such estimates are sensitive to the location of the strain gauges relative 
to the nodal lines of the initial buckle pattern. This, together with 

variations in user judgment, can lead to ambiguous answers. 

A method is presented here for determining experimental buckling 
loads for shear panels. The method employs the "diagonal tension factor", 
k, which is a measure of the amount of diagonal tension developed in the 
buckled skin. Following the model of Kuhn [1952 & 1956], which is also 
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described by Perry [1950], and referring to Figure 1, which is adapted 
from Figure 8 of Kuhn [1952], the applied shear flow q is assumed to be 

divided into two parts, a portion (1-k)q carried by shear resistant (SR) 

action of the skin, and the remainder kq carried by diagonal tension (DT) 
action in the (post)buckled skin, i.e. by a tensile stress in the direction of 
the folds of the buckled skin. The objective of the method is to derive an 
expression for k as a function of membrane stress resultants N x \ N y - and 
N x y=q>0*, where x' is normal to the stiffener direction, and y’ is parallel 
to the stiffener (or axial) direction. These resultants are calculated from 
strain gauge readings (a considerable effort in itself) and analytically 
determined skin stiffnesses. Because k is a measure of the amount of 

diagonal tension developed in the skin, a zero value of k infers initial 

buckling. Therefore, an unambiguous estimate of the experimental 
buckling load can be obtained by plotting k against load to determine the 
load at which k is zero. The expression relating k to N x -, N y - and q is 
derived next. 


From standard stress resultant transformation equations (or use of 
Mohr's circle, see Perry [1950]), the following expressions for the skin 
stress resultants, N x , N y , and N xy , in the diagonal tension axis system x, y 
(where, see Figure 1, x is parallel to the direction of the diagonal tension 
fold, and y in normal to it in the plane of the skin) can be written in terms 
of the applied skin shear flow q: 



I (l -k) q sin(2a) 1 
|-(l -k) q sin(2a) 

\ (l - k) q cos(2a) 


SR 


+ 


2kq/sin(2a)| 

0 ( 

o 'dt 


J 

1 


(l - k) q sin(2a) + 2kq/sin(2a)| 
- (l -k) q sin(2a) j 

(l -k) q cos(2a) ) 


(2), (3), (4) 


In these equations, the diagonal tension angle, oc, is the angle the diagonal 
folds make relative to the normal to the stiffener direction (i.e., a i s 
measured from the structural axis, x', to the diagonal tension fold axis, x). 
Equations (2) and (3) are given in Kuhn [1952 & 1956] and Perry [1950], and 


• Throughout this work, q is taken to be positive, as is its buckling value, q C r- 
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Eq. (4) stems directly from the stress resultant transformation equations. 
From Eq. (4) it may be concluded that the direction of the diagonal tension 
fold is generally not the same as the principal stress direction unless all 
the shear is carried by diagonal tension (k=1), or if a ■ 45°. 

To relate k to the stress resultants in the structural axis system, 
we invoke the first stress invariant, namely, that the sum of normal 
stress resultants is invariant with respect to an orthogonal 

transformation of coordinates, i.e., 

Nx + Ny = Nx' + Ny' (5) 

and we substitute N x and N y from Eqs. (2) and (3) into this equation to 
obtain 


k JN X +Ny) gin(2tt ) (6) 

2q 

To express k as a function of the three stress resultants only, we need a 
relationship between a and the resultants. To obtain this relationship, we 
equate the maximum principal stress resultant expressed in the structural 
(x\ y') and diagonal tension (x, y) axis systems to get 

(N«' + Ny) | J ] - Nyf ~ q2 (Nx + Ny) [ ^ (Nx - Nyf ~ ~T 

Eliminating the sum of the resultants on each side of this equation by 
virtue of Eq. (5,) and then squaring the radical gives 

(Nx' - Ny 1 ) 2 + 2 _ (Nx - Ny ) 2 + N 2 y 
4 4 

which is also recognizes as being the square of the radius in the x and x' 
directions of Mohr’s circle for stress resultants. To eliminate the stress 
resultants N x , N y and N xy in the diagonal tension axis system, we 
substitute Eq. (1), (2) and (3) into this equation to get 


(N X ' - Nyf 


_2kq_ + 
sin(2a) 


12 


2(1 -k)q sin(2a) 


+ 4(l -k^q 2 cos 2 (2a) - 4q 2 
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2kq ^ 
_sin(2a)_ 

2kq | 

E 

_sin(2a)_ 


+ 8k(l -k)q 2 + 4(l -k^q 2 - 4q 2 


4k q 


2 n2 


( 7 ) 


q on the right hand side of this equation may be expressed in terms of 
N X ’ and Ny- by summing Eq. (1) and (2), and using Eq. (5) to give 


(Nx- - Nyf - 


2kq 


sin(2a) 


(8) 


so that Eq. (7) becomes 

(Nx 1 ■ N y f = (N X ' + Nyf - (N x - + Nyf sin^(2cx) = (N X ‘ + Nyf cos 2 (2oc) 
Thus, the angle of diagonal tension is 

a = — cos' 1 [ Nx ' Ny | (9) 

2 \N X ' + Ny-/ 


To get the final expression for k in terms of the stress resultants in the 
structural axis system, we square Eq. (6) for k to get 

k 2 , (Nx' + Nyf sin 2 2a _ (Nx' + Nyf ( 1 . cos 2 2a ) 

4q 2 4q 2 

and substitute a from Eq. (9) to yield 


|<2 _ (N X ' + Ny') 2 (4N X ’Ny') 

4q 2 (N x - + Nyf 

from which, finally 

k _ VNx' N y - (10) 

q 
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Experimental values of k, as a function of load level, are obtained 
from this equation by calculating the stress resultants N^, N y - and q from 
strain gauge readings and analytically determined skin stiffnesses. As 
mentioned earlier, the parameter k is a measure of the amount of diagonal 
tension developed in the skin. Therefore, by plotting k against the load P 
and extrapolating the plot back to the load at which k is zero (which 
infers buckling), an unambiguous estimate of the buckling load, P cr , can 
be obtained. The corresponding value of the buckling shear flow, q cr , i s 
obtained from a (very nearly linear) plot of the skin shear flow, q (also 
calculated from the measured strains), against P, and determining the 
value of q corresponding to P cr . 

B.3 Buckling Results Based on the “Diagonal Tension Factor Method” 

Panel A1 was instrumented with back-to-back rosette gauges on the 
skin at the center of the center bay between the two stiffeners, and at the 
quarter points of the center and one side bay, see Visconti [1988]. The 
measured strains recorded during loading to failure were averaged for 
each pair of back-to-back gauges. Using calculated skin stiffnesses, the 
stress resultants N x *, N y - and q were obtained, and these were then 
averaged for the three locations The factor k was then calculated from 
Eq. (10) and plotted against load, as shown in Figure 3. By linear 
extrapolation back to the load at which k is zero, the estimated buckling 
load Per and the corresponding value of q C r are 

* LI- P C r = 37 5 K 'P S > Per = 634 PPi (11) 

The same procedures were applied to panel B1, as shown in Figure 3, 
resulting in the following values: 

• BJ.: P cr = 17.6 Kips, q cr = 302 ppi (12) 

Panel Cl was instrumented with back-to-back rosette gauges at the 
center of each bay of this one-stiffener panel. The recorded strains were 
averaged for each pair of back-to-back gauges. The value of the diagonal 
tension factor k was then calculated. This procedure was performed for 
three test runs: a preliminary run to 16 Kips in increments of 2 Kips; a 
second run to 10 Kips in 1 Kip increments to more accurately record the 
buckling behavior, and the final run in 4 Kips increments to failure (see 
Visconti [1988]). The plots of k versus load for each test are shown in 
Figure 4, together with a linear curve fit to the combined set. 
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The results for panel Cl are: 

• Ql: (P C r)k=o = 3-96 Kips, (q cr ) k=0 = 60.8 ppi (13) 

The subscript “k=0” is appended to these results for Cl to distinguish 
them from the STAGS results given in Eq. (1). Comparison of Eqs. (13) and 
(1) for Cl shows that the method developed here gives a value of q C r that 
is within 4% of the STAGS prediction. Either of the two close value of q C r 
for Cl can be used as input to SNAPPS with inconsequential differences in 
results. We choose to use the STAGS value (Eq. (1)). 

B.4 Closing Comment 

Although the above favorable correlation between values of q C r 
obtained by the “diagonal tension factor method” and STAGS i s 
encouraging for panel Cl, it is to be emphasized that more work is 
required to verify the adequacy of this new method. 


B-6 



(in direction of diagonal tension fold) 



FIGURE 1 DIAGONAL TENSION STRESS RESULTANTS (ADAPTED FROM KUHN, [1 952]). 

THE APPLIED SHEAR FLOW q IS DIVIDED INTO TWO PARTS, A PORTION (l-k)q CARRIED BY SHEAR 
RESISTANT ACTION, WHERE k IS THE "DIAGONAL TENSION FACTOR", AND THE REMAINDER kq IS CARRIED 
BY DIAGONAL TENSION ACTION. THIS FIGURE SHOWS THE STRESS RESULTANTANTS DUE TO EACH 
PART. 





0.25 



B-9 




TEST, Failure Run 



to , 't CO CM T-; O 

o o o o o o 


M ‘UOlOVd N0ISN31 TVNOOVIQ 


B-1 0 


OF THE BUCKLING LOAD FOR PANEL Cl 





APPENDIX C 


PRINCIPAL STRESS RESULTANTS IN SKIN 
C.1 Introduction 

This appendix presents expressions for the principal inplane 
membrane stress resultants Ni and N 2 as a function of the shear load. Ni 
is positive when it is in tension, and N 2 is taken to be positive when it i s 
in compression. Expressions are also given for the inplane stress 
resultants N X ’ and Ny, which are, respectively, normal to the stiffener 
direction, which is the global x' direction, and parallel to the stiffener 
direction, which is the global y' direction. The global (x\ y') axis system 
is also referred to as the structural axis system. The skin is loaded by an 
applied inplane shear flow q relative to the structural axis system. The 
shear flow can be applied in either of two directions for a fixed stacking 
sequence. Throughout this work, q is taken to be positive, as is its 
buckling value, q cr . 

Prior to buckling, the principal stress resultants equal q and make 
an angle of 45° relative to the global axes. The objective here is to 
determine their values after buckling. The expressions for the stress 
resultants are based on two simplifying assumptions (Eqs. (1) and (5)) 
that are reasonable for the objective of developing a simplified 
postbuckling method. The validity of these assumptions is confirmed by 
comparing predictions based on them with test results for the three test 
panels A1, B1, Cl described in Appendix A. 

C.2 Principal Stress Resultants 

The angle of diagonal tension (see Figure 1 of Appendix B) for flat 
shear panels is assumed to be 


a = 45°. 


( 1 ) 


Figure 1 of Section 3 shows that this is a good approximation for the 
three tested. The experimental values of (X in this figure are based on Eq. 
(9) of Appendix B, which is repeated below for completeness: 


a = — cos* 1 
2 


( Nx' - Ny' \ 
W + Ny/ 
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For this assumed value of a, Eq. (4) of Appendix B states that the 
shear stress resultant, N xy , in the diagonal tension axis system vanishes. 
Hence, the direct stress resultants, N x and N y , are in fact principal 
resultants, Ni and N 2 , i.e., the principal stress resultants line up with the 
diagonal tension folds at all load levels. Noting that N 2 is assumed to be 
positive when it is in compression, Eqs. (2) and (3) of Appendix B give 

Ni=(1+k)q (2) 

and 

N 2 = (1-k)q (3) 

where k is the diagonal tension factor described in Appendix B. Based on 
test results for aluminum panels, Kuhn [1952] gives the following 
empirical expression for k: 

k|<uhn-DT = Tanh[0.5log(q/q cr )] (4) 

For the three panels considered in the present work, Figures 2 through 4 of 
Section 3 compares principal stress resultants obtained from Kuhn’s 
diagonal tension (DT) model based on k from Eq. (4) with results from 
tests and STAGS (for Cl only). The STAGS results agree well with those 
from test, but the Kuhn-DT model slightly under-predicts the maximum 
principal diagonal tension stress resultant, Ni, but considerably over- 
predicts the minimum principal stress resultant, N 2 , in the buckled panels. 
Furthermore, as may be observed from Figure 1, the Kuhn-DT model 
predicts a significant degradation in shear stiffness with increasing load 
in comparison with the measured and STAGS results for Cl. For example, 
at the test failure load, the shear strain predicted by the Kuhn-DT model 
is about 30% higher than that of the test and STAGS for Cl. Also, 
examination of the STAGS and test results in the figure reveals that the 
overall shear stiffness of the panel is constant*, and is equal to the initial 
shear stiffness (c.f., the linear and nonlinear STAGS analysis results.), a 
result that we use in Appendix E, where we calculate the maximum out- 
of-plane displacement as a function of load. In view of these findings, we 
decided not use Kuhn’s empirical expression for k. Besides, we want the 


' Note, however, that there is a reduction in overall stiffness at the bifurcation buckling point, 
of course, as is evidenced from a plot of the overall corner displacement, S, at the loaded 
comer against the fixture load, P, as shown in Figure 7 of Sobel and Sharp [1994]. The slope 
of the P-S curve is a measure of the overall stiffness of the panel and therefore represents all 
stiffnesses, including out-of-plane bending stiffnesses, as is evident from the energy 
consideration that PS/2 equals the panel’s total strain energy (bending plus membrane). 
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SNAPPS methodology to be devoid of all empiricism. To obtain k, and 
hence Ni and N 2 via Eqs. (2) and (3), we make a further assumption; 
namely, that the compressive principal resultant after buckling, N 2 , is 
constant for all values of q and equal to the buckling shear flow, q cr , i. e., 

N2 - q cr > ( for Q * Per) ( 5 ) 

With this assumption, Eq. (3) gives 

k = 1 - q cr /q (6) 

Inserting Eq. (6) into Eq. (2) provides the following expression for the 
tensile principal resultant after buckling: 


N 1 = 2q-q cr , (for q > q cr ) (7) 

This view of events is similar to that proposed by Wagner in the "frame 
analogy", as described by Kuhn [1952] and illustrated in Figures 3a and 3b 
of Kuhn. It is noted that Kuhn first gives an expression for k (Kuhn’s Eq. 
(26)) that is the same as the “Wagner-type” expression for k ((Eq. (6) 
above), but later abandoned it in favor of an empirically determined 
diagonal tension factor, Eq. (4). 

It is important to point out that if N 2 was assumed to increase (in 
an absolute sense) above its value at buckling, then Ni would have to 
decrease by the same amount in order for the sum of (N-|+ N 2 ) to remain 
constant (in accord with the first stress invariant principle), and this 
would lead to a higher predicted failure load, qfail- Thus, the assumption 
that N 2 = q cr after buckling is conservative. On the other hand, it would be 
unrealistic to assume that N 2 is zero for all values of q, i.e., to assume a 
state of pure diagonal tension. 

C.3 Direct Stress Resultants In Structural Axis System 

Using standard equations for transformation of stresses from one 
orthogonal axis system to another, we readily obtain from the expressions 
for the stress resultants in the diagonal tension axis system (Eqs. (5) and 
(7)), the following expressions for the stress resultants in the structural 
axis system: 

N X ' = Ny' = q-q cr , (for q > q cr ) (8) 
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As a check we substitute these equations into the general expression for 
k, (Eq. (10 )) of Appendix B, repeated below for convenience, 


k _ VNx' Ny' 

q 

to obtain 

k = 1 - q cr /q 


(9) 


which is Eq. (6). 


C.4 Closing Comments 

The principal stress resultants described by Eqs. (5) and (7) are 
averages over the extent of the buckle. However, they can be used to 
approximate integrated quantities such as the buckle wavewidth, c 
(Appendix D); the maximum out-of-plane displacement, Smax (Appendix E); 
and the rotation at the end of the tie-rod, 0 (Appendix F). 

Postbuckling values of the principal stress resultants are also 
needed for local considerations at the toe of the flange. There the skin i s 
flat and can be assumed to be in a state of pure shear. This means that 
the stress resultants in the skin under the toe in the diagonal tension 
field direction are given by 


Nnoe = N 2 toe = q. (for all q) (10) 

It is to be recalled that N 2 toe is taken to be positive when it is in 
comression, and that q>0. Equation (10) is used in Appendices H, I, and K 
to describe local stresses in the unbuckled skin under the toe of the 
flange. 
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APPENDIX D 


BUCKLE WIDTH & LENGTH 


D.1 Buckle Width 

Figure 1 depicts a buckle bounded by adjacent nodal lines and the 
stiffener flanges. The buckle width c is determined next from 
considerations of out-of-plane equilibrium. If we assume that the buckle 
may be approximated by the rectangle ABCD of length l and width c, then 

we may take the out-of-plane displacement 5 as 

5(^,ri)=5 m axsin^-sin^- ^ 

where the axis system has its origin at B; \ is measured from B to C, rj i s 
measured from B to A. The stress resultant in the ^direction is tensile 
and is denoted by Ni ; that in the ri direction is compressive and is denoted 
by N 2 (a positive N 2 is compressive). As discussed in Appendix C, they are 
assumed to be constant throughout the area bounded by the rectangle 
ABCD. This piece of the buckled sheet has to be put into equilibrium 
normal to its original flat plane, under the components of N-| and N 2 acting 
through the slopes along the edges of the distorted shape. 

The total out-of-plane force due to the tensile stress resultant Ni 
acting at the angle d8/d% along one edge of the buckle (such as AB in 
Figure 1) is then given by 



dr| = 2 Ni 


$max ^ 


This must be balanced by the following out-of-plane force due to the 
compressive stress resultant N 2 acting at the angle d5/drj along a nodal 
edge of the buckle (such as BC in Figure 1): 

ci 

N 2 W d£, = 2 N 2 5 m ax £ 

\dn/Ti=0 c 

JO 
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By equating these two forces, we obtain 

2 Ni 8max ^ = 2 N 2 5max ^ 
y 0 

It is seen that 5 max in this equation conveniently cancels out (8 max is 
determined in the next appendix) thereby providing the following 
expression for the buckle aspect ratio: 



In terms of known quantities, and assuming that the diagonal tension 
angle a is 45°, the length i can be written as 

t = Mk) - c 
\coscc / 

From Figure 1, is the length OC, denoted by L, the known toe-to-toe 

cosa 

distance between flanges measured in the diagonal tension axis system, 
i.e., 



(3) 


Thus 


i = L - c 

so that Eq. (2) yields the following expression for the buckle wavewidth: 



in which (see Appendix C) 

Ni=2q-q cr , (q>q cr ). 
N2= Per. (q>Pcr)- 
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It is to be recalled that positive Ni is tensile, and positive N 2 is 
compressive. 

At buckling, for which N-|= N 2 , Eq. (4) yields 



whereas, from Timoshenko and Gere’s [1961, page 383] buckling analysis 
of a long, simply supported, isotropic plate under shear, it may be shown 
that 


(?) =5-^ 

'L./cr.Timo ^ 

It is proposed that Eq. (4) be modified to pass through this buckling point 
to give 



(7) 

From the last term on the right, it is seen that the wavewidth c is 
expressed as a function of the load level q. Effectively, we determined 
the shape of the c vs. q curve from equilibrium considerations, and 
adjusted the amplitude (by 22.5%) of the curve to have the curve pass 
through the buckling point. 

D.2 Buckle Length 

The axial wavelength, Lb, which is the projected length of the 
wavewidth along the stiffener axis, is used in Section 3 for correlating 
SNAPPS results with those from test and STAGS. From Figure 1, it is seen 
to be given by 


Lb = 


e 

cos a 


( 8 ) 


From Eqs. (3) and (7), and with the angle of diagonal tension equal to 45°, 
the wavelength may be written in the following dimensionless form: 
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Nfc 

Ni 


1 + 


N2. 

Ni 


(9) 


-^ = 2 VTT 5 

bn 



2.45 



where 

b tt = b-2L f (10) 

is the toe-to-toe width of the skin, measured normal to the stiffener 
direction (see Figure 1). Note that this form for the wavelength is valid 
for all panel geometries. A comparison of predicted and measured 

wavelengths based on Eq. (9) is given in Section 3. 
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APPENDIX E 


MAXIMUM OUT-OF-PLANE DISPLACEMENT 

As we saw in the preceding appendix, the out-of-plane displacement, 
5, through its effect on the slopes at the boundaries of the buckle, governs 
the out-of-plane equilibrium of the buckle, from which we were able to 
determine the buckle width c. Now we need to find an expression for the 
maximum value of the displacement, 5max- The following derivation is 
based on the kinematic argument that the buckled sheet must fit into the 
distorted shape imposed by the nodal lines and the stiffeners, and that, 
after buckling, it does this by moving out of its original flat plane into the 
characteristic shear buckling pattern. 

Figure 1 shows a single buckle, bounded by adjacent nodal lines and 
by the toes of the stiffener flanges. The length of the nodal line OC i s 
assumed to remain unchanged in a shear distortion, and the length d of the 
shorter diagonal AC is given by 

dj 2 = L 2 + Lb 2 - 2LLbCOs(7r/2-<J>) 

= L 2 + Lb 2 - 2LLbsin<t> for small values of y. 

Imposing a shear strain y, such that <|> becomes <t>+y, the new length of the 
diagonal AC is given by 

d n 2 = L 2 + Lb 2 - 2LLbsin(<j) +y) 

= (L 2 + Lb 2 - 2LLb sin(>) - (2LLbY cos<J>) for small values of y. 

Increasing the shear strain from the buckling value y C r to y, and expressing 
the change in the length of the diagonal AC as a strain £d [i.e. where e<j = 
(new length - old length)/old length] gives 


£d = (dn _ dj)/dj 


L 2 +Lb 2 -2LLb sin({j 

]1/2 

i -2LLbycos(j)J 

r o o i 1/2 

L^+Lb -2LLbSin(J) -2LLbY C rCOS<}) 


L 2 +Lb 2 -2LLbSin(j) 

-2LLbY cr cos(]) 

1 / 2 


By expanding the radicals, and neglecting terms in y 2 or higher powers, the 
expression for £d can be approximated by 
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- LLb(Y-Ycr) gos<j> 

L 2 + Lb 2 - 2LLb sintj) 

But L= b n sec(|) (see Figure 1) and = (L^ + Lb^ - 2LLb sin<})) so that the 
approximation for 8<j may be written as 


Ed = - ^(Y-Ycr) 
d 2 


( 1 ) 


where the minus sign denotes compression. 

Now, to obtain an alternate expression for ed in the postbuckled 
panel, define an axis X along the diagonal AC, and assume the out-of-plane 
displacement to be 



It is straightforward to show that this assumed displacement variation 
along the diagonal d is consistent with the double-sine displacement 
shape (Eq. (1) of Appendix D) used in the preceding analysis for the 
determination of the wavewidth c. Later we will show that 5 ma x is 
insensitive to the assumed displacement shape. Integrating the length ds 
along the curve from 0 < X < d, using 


-dS- = 1 + =1+1 j rc^max | 2 s j n 2|27cx| f or sma || 

dx v Ux) 2 1 d I ' d ' < 

gives the length the total S measured along the curve as 


max 


d 1 + 1 
\ 4 


ft Smax l 2 


d I 


Expressing the change of length as a strain e<j [= (d-S)/d ] yields 

1 /jtSr ' 2 


d - d 1 


^d = 


+ 1 | rc$max 
4 I d 


'max 

2d 


( 2 ) 
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Equating the two expressions for ed, Eqs. (1) and (2), provides the 
following closed-form expression for the maximum out-of-plane 
displacement: 

§max = ^Vl_t>btt(Y-Yc r) (3) 

It may be noted that the precise shape of the out-of-plane distortion 
is not very critical, because we are integrating along a curve and 
expressing the result as an average (ed)- In a study in which different 
expressions for the assumed shape were used, such as circular arcs, 
parabolas, fourth-order polynomials, sine and sine squared waves, Eq. (3) 
retained the same form and the multiplier in front of the radical changed 
by less than 7% from the higher and hence more conservative value of 2 In . 

Because the STAGS analysis and test results for panel Cl [Sobel & 
Sharp, 1994] and the test results for panels A1 and B1 [Visconti, 1988] 
revealed that the skin shear stiffness Gt did not change with applied load 
level (see Figure 1 of Appendix C, and discussion therein), Eq. (3) can be 
written as 

6 max= lV Lt>b "( 9 &rf <4) 

In this expression, b tt is the toe-to-toe distance measured normal to the 
stiffener direction, i. e., 

b tt = b-2Lf (5) 

where b is the stiffener spacing, Lf is width of the flange measured from 
stiffener centerline to toe, and Lb is the buckle wavelength projected 
along the stiffener axis. An expression for Lb is given in the preceding 
Appendix. 

An alternative form for the maximum out-of-plane displacement 
may be obtained by inserting 

Lb = c sec(|) and b„ = L cos(j) (see Figure 1), 

into Eq. (4). This yields the following formula, written in terms of 
diagonal tension axes: 

5max = (6) 
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In this equation, L is the toe-to-toe length of the diagonal tension fold, 
and c is the buckle width for which the following expression is derived in 
the preceding Appendix: 



(7) 


The form for the maximum out-of-plane displacement given by Eq. (6) is 
used in the SNAPPS code of Appendix N. 

From the last two equations, the maximum out-of-plane 
displacement may be written in the following dimesionless form: 




This dimensionless form for 5max is used in Section 3 to compare 
predicted and measured values of the maximum out-of-plane 
displacement. Note that form is valid for all panel geometries 
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APPENDIX F 


TIE-ROD ANALYSIS 


F.1 Introduction 

The development of the simplified analysis now proceeds to the 
"tie-rod" analysis in which a unit width of skin parallel to the long axis of 
the diagonal tension buckle is assumed to be loaded axially by the tensile 
stress resultant, Ni, and by a lateral load p, which is normal to the 

original unbuckled plane of the skin and which represents the out-of-plane 

component of the compressive stress resultant N 2 . For this simplified 
analysis, it suffices to assume that Ni and N 2 are constant along the 

length of the tie-rod, but, in fact, both stress resultants vary along the 

length and fall to the unbuckled values at the toe of the flange. Because 

the angle through which N 2 acts varies, the lateral load p is assumed to 

vary sinusoidally along the length of the nodal line, a variation consistent 
with those used in the wavewidth (Appendix D) and maximum out-of-plane 
displacement (Appendix E) analyses. The objective of the tie-rod analysis 
is to obtain a closed-form expression for the rotation 0 at the end of the 
tie-rod in terms of the axial load Ni, the maximum out-of-plane 

displacement 8 ma x. and the moment M at the end of the tie-rod. 

Figure 1 shows the tie-rod model. The upper half of Figure la 
depicts the assumed bi-sinusoidal buckle pattern in a typical diagonal 
tension fold, and the shaded strip in the lower half represents one-half of 
the tie-rod model, which extends along the entire length, L, of the buckle*. 
The tie-rod model is further displayed in Figures 1b to Id. From Figure 
1c it is seen that the vertical (out-of-plane) component of the 

compressive membrane force, N 2 , produces transverse shear forces at the 

edges of the tie-rod. As just mentioned, these edge forces are assumed to 
vary sinusoidally in the direction of the fold (Figures 1b and Id) and are 
represented by a lateral (out-of-plane) running load, p (Figures 1b and 1c). 


* The stiffener is not shown in the figure, to avoid excessive clutter. However, the downward 
deformation of the buckle, implies that the stiffener is on the top surface of the skin (on the 
surface with the negative y coordinate. 
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The running load p is reacted by the edge forces R (Figure Id). Thus, R 
balances the vertical component of N2 that varies along the longitudinal 

edges of the strip. 

For simplicity, the tensile axial load in the tie-rod model is held 
constant at it's average value (i.e. Ni = 2q-qcr. see Appendix C) and the 
solution is then tractable. In fact, the axial load must be slightly above 
the average value over most of the length of the tie-rod, but must fall, 
over an undefined— but assumed small — distance at each end, to the value in 
the unbuckled sheet (i.e. Nitoe = see Appendix C ) 


F.2 Solution Of Tie-Rod Problem 

Figure 2 displays the notation and sign conventions. The x axis is in 
the direction of the diagonal tension fold, and it's origin is at mid-length 
of the fold. The y axis is normal to the plane of the skin. Because the 
solution of the tie-rod problem is straightforward, it is presented in 
more-or-less outline form. 

. Relate Reaction Rto Magnitude (p n ) of the Pressu re P 



• Differential Equation (DEV [Timoshenko. 19561 

D sr w - N-| w" = p 0 cosj^-J ( 2 ) 

where ( ) = -^-(x) and D$r is the bending stiffness of the skin (subscript s 
dx 

for "skin”, subscript "r" for "rotated" in the direction of the diagonal 
tension fold). 

Letting: 

x 2 = Ni (3) 

Dsr 
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The DE becomes: 


w’ ,H - X 2 w" = -£°-cosM (4) 

D sr ' L / 

• Complementary Solution. y j c 

The solution of the homogenous form of Eq. (2) is 
w c = ACosh (Xx) + BSinh (\x) + C +Dx 

• Particular Solution. w P 

A particular solution of Eq. (2) is easily shown to be: 

w p = r cos j^-J (5) 


where r, which has dimensions of length, is given by 


r = 


L 

_ZL 


1 + |-*-f 

XU 


(6) 


• Total Solution, w = w c + w p 

w = ACosh (?ix) + BSinh (>.x)+ C +Dx + r Jp- cosj^-j (7) 

Because the displacement must be symmetric about the midlength, x = 0, 
the two asymmetrical terms (B Sinh (Xx) and Dx) must vanish. 

Therefore, 

w = ACosh (Xx) + C + r Jp- cos (8) 
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The arbitrary constants (A and C) in this equation are determined from the 
following boundary conditions (BC) at x = ±^): 


BC1: M(±^) = M, 


w (± — ) = 0 

2> 


(where M is depicted in Figure Id) 


Using M = -D sr w" (see Figure 2b) and Eq. (3), BC1 becomes: 
w" l ± L) = -M_ = -x 2 M. 


D 


sr 


Ni 


Substituting Eq. (8) for w into the above expression yields the arbitrary 
constant A: 


A = -SechW-^- 
V 2 / Ni 


(9) 


The constant C is obtained from BC2 and Eq. (9): 

C = -ACosh (XL) = +Sech W Cosh = -M- 

\ 2 I 12 / \ 2 / Ni Ni 


( 10 ) 


Substitution of these constants into Eq. (8) provides the solution for w: 


w(x) = 


1 - 


j ! ' — 

• Rotation at Left End x = - — , 0 

t 2 / 


Cosh (Xx) 

Cosh j 

N 


12 /J 

t End 

(*- 


_M_ + r _B_ cos (— 
Ni Ni \ L 


(11) 


Differentiating Eq. (11) once gives: 

s- sin (ax.) 


w - (x) .. x iiSinhM. r-B-a 


Ni Cosh W Nl 
2 / 


( 12 ) 
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• Tanhf^H term in 0 
2 I 


Tanh W -> 1 for > 2, or XL > 4 

’ 2 1 2 

• 1 - Sechj^-J term in 5 ma x 


Assuming a 10% error in neglecting Sech|^J in comparison 
with unity gives 


Sech W = 0.1 => Cosh W = 10 

1 2 / 12 / 


=>^-=3=»XL>6 

2 


• 1 + (-S-) 2 term in denominator of r (Eq. (6)) 
\XU 


Assuming a 10% error in neglecting MM 2 in comparison with 

'XL/ 

unity gives 


(JL.) 2 = o.l => XL = -5= = 9.9, or XL > 10 

'XL/ VoT 
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j 1 1 ii 


=>r = L 

Ti 

With the above approximations, the equations for the rotation and 
maximum out-of-displacement simplify to 



(18) 

Ni Ni 


8max 88 

(19) 

Ni * Ni 



In view of the above discussion the simplified equations for 0 and 8 m ax are 
valid for: 

XL >10 ( 2 °) 

The simplified expressions given by Eqs. (18) and (19) are used in this 
report. 
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FORCE, N2, PRODUCES TRANSVERSE SHEAR FORCES AT THE EDGES OF THE TIE-ROD. THESE EDGE FORCES 

ARE ASSUMED TO VARY COSINUSOIDALLY IN THE DIRECTION OF THE FOLD (FIGURES b AND d) AND ARE 
REPRESENTED BY A LATERAL (OUT-OF-PLANE) RUNNING LOAD, p (FIGURES b & c). THE RUNNING LOAD p IS 
REACTED BY THE EDGE FORCES R (FIGURE d). THUS, R BALANCES THE VERTICAL COMPONENT OF N 2 THAT 

VARIES ALONG THE LONGITUDINAL EDGES OF THE STRIP. 


in direction of diagonal H UA in +v dir on positive face (+x coord) 
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APPENDIX G 


SKIN BENDING MOMENT AT TOE, M toe 


G.1 Introduction 

The purpose of this appendix is to determine the bending moment in 
the skin immediately under the toe. This is done by linking the "tie-rod" 
model, described in the previous Appendix F, to a model of the combined 
stack of the skin, flange and adhesive layers between the toe of the flange 
and the heel (or stiffener centerline). The latter model is based on that of 
Tsai [1983] in that it is cut from the stiffener in the same diagonal axis 
as the tie-rod model of the skin, but departs from Tsai’s model in that i t 
includes the entire stack of layers, adhesive and SynCore of the skin and 
the flange (for this reason, properties of the total stack have the 
subscript "tot" and, because they are in the diagonal tension or rotated 
axis system, "totr"). The model of the combined skin and flange is 
represented by a cantilever in the diagonal tension direction, clamped at 
the stiffener centerline and loaded at its tip by the transverse shear force 
and bending moment from the end of the tie-rod. The moment is 

determined by requiring rotational compatibility between the tie-rod end 
and the tip of the cantilever. Because the transverse shear force i s 
related to the moment through the maximum out-of-plane displacement 
Smax of the tie-rod, the only unknown in the equations for the rotations is 
the moment. 


G.2 Rotation at End of Tie-rod (see Figure 1 of Appendix F) 

From Eq. (18) of Appendix F, the rotation at the left end of the tie- 
rod, is given by 


0 = X_M + _B_ 

Nt Ni 

and from Eq. (19) of the same appendix, the transverse shear force at the 
end of the tie-rod is 


( 1 ) 


_EL - 2L 
Ni L 


•'max 


M.) 

Nil 


In these equations, 


(2) 
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(from Eq. (3) of Appendix F) (3) 



L is the toe-to-toe length of the buckle, see Figure 1 of Appendix E, 
5max is the out-of-plane displacement at the center of the tie-rod, 
as described in Appendix E, 

Ni is the membrane stress resultant in the buckled skin in the 
diagonal tension direction, 

D S r is the bending stiffness of the skin about it’s neutral axis 
normal to the long axis of the buckle. 

Eliminating R/Ni from these equations gives 

9 = (4) 

Ni L \ Ni / 

G.3 Rotation of "Flange and Skin" Cantilever Model 

Figure 1 shows the beam model of the combined skin and flange. The 
beam is assumed to be clamped at the heel, an assumption justified by the 
STAGS results of Figure 2. Any elementary strength of materials text 
provides 


Qtotr = 


1 + (M, oe - qz) 

2 Dtotr Dtotr 


(5) 


where 

Wtoe = the transverse shear at the tip of the cantilever, 

Mtoe = the bending moment at the tip of the cantilever, 

Dtotr s the bending stiffness of the combined stack of skin, adhesive 
and flange, 

Lfr = the heel-to-toe width of the flange measured in the direction 
of the diagonal tension fold, 

qz = the moment due to the offset of the skin stress resultant from 
the midplane of the combined stack, as described below. 

As described in Appendix F, the tensile axial load in the tie-rod 
model is held constant at it’s average value (i.e. Ni = 2q-q cr , see Appendix 
C) . In fact, the axial load must be slightly above the average value over 
most of the length of the tie-rod, but must fall, over an undefined — but 
assumed small — distance at each end, to the value in the un bu ckl e d sheet 
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at the toe (i.e. Nitoe - q, see Appendix C). When the tie-rod is connected to 
the tip of the cantilever the change in Ni at each end of the tie-rod must 
betaken into account. The tensile end load transferred from the end of 
the tie-rod to the cantilever is therefore Nn oe = q and njal Ni =2q-q C r- At 
first sight, it appears that equilibrium is being violated, but the change is 
assumed to occur over a short length at each end of the tie-rod and is 
brought about by local in-plane shear stresses in the skin near the toe. 
These local changes are assumed not to invalidate the solution of the tie- 
rod problem. 

In Eq. (5) an additional moment (qz) has been added to Mtoe to 
represent the moment due to the membrane stress resultant in the 
unbuckled skin at the toe (i.e. Nitoe = q) being offset a distance z from the 
centroid of the combined stack of the skin and flange layers (see Figure 1). 
z is given by 


Etfits + t a +-^- + Et S pf- 

z 1 2J 2_ - k_ (6) 

Etfr + Et sr 2 

where 

Etf r = membrane stiffness of the flange in the diagonal tension 
direction 

Et S r = membrane stiffness of the skin in the diagonal tension 
direction 

tf = flange thickness 
t s = skin thickness 
t a = adhesive thickness 

G.4 Determination of Bending Moment & Transverse Shear Force at Toe 

By enforcing the following conditions, obtained by comparing Figure 
1 of this Appendix to Figure Id of Appendix F: 


M = - Mtoe 

(7) 

R = + W t0 e 

(8) 

0 = + 0totr 

(9) 


we obtain the following expression for the bending moment in the skin at 
the toe: 
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7C 5max | -j 

Mice =7-^ 

X. - rc/L 
l Ni 


Lfr Ni | | Lfr q z 

2 Dtotr l Dtotr 

Lfr , 71 L fr | 

Dtotr 2 L Dtotr' 


( 10 ) 


where 




c = 

Vi. 5 

L "V 

Ni ' 

i (from 

Eq. 

(7) 

of 

Appendix 

D) 

(11) 



1+^ 

/* : 









1 ^ 

/ Ni | 








Smax - 

71 

iff 

- Per) 
Gt 1 

(from 

Eq. 

(6) 

of 

Appendix 

E) 

(12) 

Ni 

= 2q 

■ Per 


(from 

Eq. 

(7) 

of 

Appendix 

C) 

(13) 

n 2 

= qcr 



(from 

Eq- 

(5) 

of 

Appendix 

C) 

(14) 


and 

Gt « the in-plane shear stiffness of the skin referred to axes 
parallel and normal to the stiffener direction. 

Also, from Eqs. (2), (7), and (8), we obtain the following expression for the 
transverse shear force in the skin at the toe: 


Wtoe = (N-|5max + Mtoe) ("15) 

Equation (10) is the full version of the expression for Mtoe and is 
used in the code listed in Appendix M. 

G.5 Parametric Study 

In the course of applying SNAPPS, it was found in some cases that 
the moment in the skin at the toe reaches a maximum and then decreases. 
The following investigation of this behavior brought to light two 
dimensionless parameters which govern the magnitude of the maximum 
value of Mtoe and the shear flow q* at which it occurs, and which have 
practical design consequences. 
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Examination of typical magnitudes of the various terms in Eq. (10) 
reveals that Mtoe may be roughly approximated by dropping terms in the 
numerator and denominator. The expression for Mt 0 e then simplifies to 


Mtoe = 


ftSmax 

[l L <r N * ) 

L 

L 2 Dtotr/ 


X 


Ni 


(16) 


By substituting for 5 ma x. c, and Ni and X in terms of q and q cr , it can be 
shown that Mtoe can be expressed in terms of two dimensionless 
parameters and a function of q/q C r- 


AjSIT - ( 
Per V D S r V 


1 - 


Pcr l-f r 
Dtotr 



(17) 


The function f(q/q or ) rises monotonically against q/qcr, but the 

presence of the parameter SsiA makes it rea ch a maximum and then 

Dtotr 

fall, as illustrated in Figure 3 where Mtoe is plotted against q for panel Cl. 
(Neither A1 or B1 reach the predicted maximum value of Mt 0 e. as is shown 
in Figures 4 and 5.) 


The form of Eq. (17) indicates that, for a specific panel, if we plot 
Mtoe against q, and search for the maximum value of Mtoe. and the value (q*) 
of the shear flow at which it occurs, then we can plot both the maximum 

moment parameter Moe. rGl ) and the shear flow q*/qcr at which i t 

Per V D sr 

,2 

occurs, against ^£ 1 — fr 
Dtotr 

Analyses, using the full version of the expression for Mtoe given in 
Eq. (10), were done for the three test panels, A1, B1, and Cl. In each case, 

.2 

several different values of ^ ^ were used, obtained by varying the 

Dtotr 

bending stiffness of the combined skin and flange stack, Dtotr- The 

2 

resulting values of and q*/qcr are plotted against — — - in 

Per V Dsr Dtotr 

Figures 6 and 7. In each figure the solid symbols denote the nominal panel 
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and the open symbols are for the variations. Both figures show that the 
results from all three panels can be approximated by a single curve, as 
predicted by the approximate form of Eq. (16). It is apparent that both 
curves are fitted quite well by hyperbolae, and that the maximum toe 
moment can be approximated by 


Mtoe.max = 0.35 (Dtotr/Lfr^) VDsr/Gt 
occurring at a shear flow 


(18) 


q* = 0.59 Dtotr/Lfr 2 


(19) 


G.6 Design Implications 

From a design point of view it is desirable to reduce the toe moment 
as far as possible because the peel stress between the skin and the flange, 
the membrane stress in the critical skin ply, and, therefore, the maximum 
principal tensile stress depend on Mtoe- The above results show that the 
maximum value of the toe moment increases as the skin/flange bending 
stiffness Dtotr is increased, or as the flange length Lf r is made shorter. 
The term VD S r/Gt is approximately proportional to the skin thickness and 
consequently the toe moment also increases as the skin thickness ts i s 
increased. 

In summary, these results confirm the need for compliant flanges 
which impose little restraint on the deformation of the buckled sheet 
(however, the stiffener must be sturdy enough to prevent the diagonal 
tension folds from progressing across the stiffener centerlines). 
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FIGURE 2 STAGS RESULTS FOR THE OUT-OF-PLANE DISPLACEMENT 
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FIGURE 4 VARIATION OF SKIN BENDING MOMENT AT TOE 
OF FLANGE WITH APPLIED SHEAR LOAD. PANEL A1. 
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FIGURE 7 SHEAR FLOW, q*, AT WHICH THE MOMENT IN THE 
SKIN AT THE TOE, M , IS A MAXIMUM 






APPENDIX H 


INTERLAMINAR SHEAR STRESS BETWEEN SKIN AND FLANGE 

In this appendix, an expression is presented for the shear stress 
distribution along the interface between the skin and the flange, based on 
the "shear lag" analysis of Kuhn [1956]. The distribution is used in the 
peel stress analysis, described in Appendix I, and the maximum value of 
the shear stress, which occurs under the toe of the flange, is used in the 
failure criterion in Appendix L. 

Figure la shows a section of the combined skin and flange model in 
the x-z plane, where x is in the direction of the diagonal tension fold, and 
z is normal to the plane of the skin. The section is of unit width in the 
inplane direction normal to the fold (i.e., of unit width normal to the plane 
of the paper). Following Kuhn, the skin-flange model is idealized as two 
axial load carrying bars joined by a shear carrying web, as indicated in 
Figure 1b. The boundary conditions are that at the toe, x = 0, the bar 
representing the skin is loaded by an axial tensile force P, and the bar 
representing the flange has zero load. A further assumption is made that 
the width of the flange Lfr is long enough for the shear stress to die away 

to zero at the heel. This assumption leads to Kuhn's Eq. (4.10) 
corresponding to Case 4 of Table 4.1. Adapting that equation to the Kuhn- 
type model of the layered model (Figure 1b) leads to the following 
equation for the distribution of the interlaminar shear stress between the 
skin and the flange: 

Kx)-E&( . EfA l-- )e-Kx (1) 

t \Ef Af + EsAs/ 

where (Kuhn’s Eq. (4.4)) 

K 2 = St/ — 1 — + _J — j (2) 

b \EfAf E S A $1 

In applying Eqs. (1) and (2) to a unit width w of the layered model, 
the axial stiffnesses EfAf and E S A S of the flange and skin in the Kuhn-type 
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model are replaced by wEtfr and wEt sr respectively (where Et denotes the 
product of modulus E and thickness t; the subscripts f and s denote the 
flange and skin, respectively, and the subscript r denotes that these 
quantities are measured in the [rotated] diagonal tension direction). The 
load P is replaced by wNitoe. (where Nitoe is the end load in the skin 
immediately under the toe, as described in Appendix C), and the web 
thickness t is replaced by w. Further, because G is not uniform throughout 
the model, the quantity G/b in Eq. (2) must be replaced by a "modulus of 
foundation", denoted by <})xz> which is given by Eq. (2) of Appendix J. With 
these substitutions, w nicely cancels out, and Eqs. (1) and (2) become 


and 


T(x)= Nitoe K 


f Et fr 
\Etfr + Et 



e- Kx 


K 2 


= <t>xz 


1 - 3 - + 

'Etfr 



(3) 

(4) 


The maximum value of the shear stress occurs at the toe (x = 0), and 
is simply denoted by Txz (without any subscripts “toe ), i.e., 


T X z =(T(x))x = o = Nitoe K 


( Etfr, 
lEtfr + Et 



(5) 


In the interest of keeping the model as simple as possible, a more 
exact analysis is not considered worthwhile at this time. However, the 
use of SUBLAM [Flanagan, 1993] to verify the distributions and magnitudes 
of the peel and shear stresses between the flange and the skin, should be 
part of any future work. 
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FIGURE 1 KUHN-TYPE SHEAR LAG MODEL FOR OBTAINING INTERLAMINAR SHEAR 
STRESS BETWEEN SKIN AND FLANGE 





APPENDIX I 


PEEL STRESS ANALYSIS 


1.1 Introduction 

A major contributor to stiffener disbonding is the peel stress 
between the skin and the flange immediately under the toe of the flange. 
Typical of the basic strength-of-materials approach adopted in SNAPPS, a 
simple beam model is used here to determine this stress. Figure 1 shows 
the model in which the skin and flange are idealized as beams that are 
coupled by an "elastic, Winkler-type, foundation" of modulus <j) z 

(determined in Appendix J). The peel stress is given by 

CpeelOO = <j>z[Wf(x)-w s (x)] ( •] ) 

where Wf and w s represent out-of-plane displacements of the flange and 
skin, respectively. The peel stress and displacements vary along the 
length of the coupled beams, as indicated in Eq. (1) by the axial coordinate 
x, which is measured from the toe of the flange to its heel, in the 

direction of the diagonal tension fold. The other coordinate used in this 

analysis is z, which is normal to the plane of the skin. The objective of 

the analysis is to obtain a closed-form expression for the peel stress at 
the toe of the flange. This stress is simply denoted by a z , (avoiding 
subscripts such as “toe”). That is, 


<5z - (<*peel)x=o 


The beams are assumed to be of length Lf r in the direction of the diagonal 
tension fold, and of unit width in the inplane direction normal to the fold. 
Lf r is the “rotated” flange (subscript “r” for rotated) width measured 
from the heel to the toe of the flange in the direction of the fold. The 
objective is to determine Wf and w s , and hence <|) z , as a function of the 
following three loadings shown in Figure 1: the transverse shear force 
Wfoe in the postbuckled skin; the bending moment Mfoe in the skin; and the 

interlaminar shear stress distribution t(x) at the skin-to-flange interface. 
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The interlaminar shear stress, whose presence complicates the analysis, 
arises from the diagonal tension field pull Nn 0 e unloading some of its load 
from the skin into the flange (see Figure 1). This shear stress causes a 
moment on each beam because it acts at a distance from the midplane of 
the flange and the skin. Thus, the present analysis may be regarded as 
being an extension of the classical "beam on an elastic foundation" method 
[Hetenyi, 1946] to the case in which two beams are elastically coupled 
together and are subjected to an offset shear traction on their interface 
surfaces. 

We next make a reasonable assumption that greatly simplifies the 
analysis, and which is consistent with our objective of obtaining a 
closed-form expression for the peel stress at the toe, which is where the 
peel stress peaks and where failure has been observed in our tested shear 
panels. The assumption is based on results from detailed analyses [Cacho- 
Negrete, 1978] for thin laminates joined by a thin adhesive layer. The 
detailed results reveal that the peel stresses between the skin and the 
flange are highly localized under the heel and toe. Because they die-away 
so rapidly from each end of the flange, the two stress concentrations do 
not interact and are essentially uncoupled. Thus, for purposes of 

determining the peel stress at the toe, we may assume that the length 
(rotated flange width, Lf r ) of the coupled beams is long enough so that 
conditions at one end (heel) do not affect the other end (toe). 
Analytically, it can be shown [Hetenyi, 1946] that this assumption is 

valid provided that the following equation involving the characteristic 
wave number (5 is satisfied: 

(3l_fr>7L 

This condition is satisfied by our tested shear panels. For example, for 
panel Cl, Figure 2 clearly shows that the rotated flange width 


Lf r = (0.75")/ cos(45°) = 1 .06" 

is "long enough" so that waves emanating from one end die out rapidly 
before they reach the other end. The decaying functions plotted in the 
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figure represent terms in the complementary solution for a beam on an 
elastic foundation. The assumption that conditions at one end of the 
flange do not affect conditions at the other end greatly simplifies the 
solution of the coupled beam equations and leads to an explicit closed- 
form formula for the peel stress at the toe in terms of the transverse 
shear force and bending moment in the skin, and the interlaminar shear 
stress between the skin and the flange. This formula requires an 
expression for the shear lag stress distribution x(x) as a function of 
distance along the width of the rotated flange. This expression is given in 
the preceding appendix, and is repeated here for completeness: 


x(x) = Nnoe 


K I — 


lEtfr+Etsi 


e 


Kx 


( 2 ) 


where 


Nitoe is the end load in the skin immediately under the toe. As 
described in Appendix C, Nitoe = d- 


K = J 4>xz I- 1 - + -M 

V lEtfr Ets,/ 

and 

(|>xz is the foundation modulus corresponding to the interlaminar 
shear stress, a closed-form expression for which is given in the 
next appendix, 

Et S r = membrane stiffness of the skin, 

Etfr - membrane stiffness of the flange. 

Et S r and Etf r are measured in the direction of the diagonal tension fold. 


1.2 Derivation of Equation for Peel Stress 


Figures 3 and 4 display the notation, sign conventions, differential 
equations, and boundary conditions for the coupled beam analysis. The 
differential equations and boundary conditions differ from those of 
Hetenyi [1946] by the inclusion here of the interlaminar shear stress, x, 
between the skin and the flange. They were derived here by writing 
equilibrium for the general differential element shown in Figure 3a, and 
using the constitutive equations given in Figure 3b. Specific forms of the 
differential equations and boundary conditions for the skin and flange are 
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summarized in Figure 4. It to be observed from that figure that the 
interlaminar shear stress appears in the boundary conditions, as well as 
in the differential equations, and that the differential equations are 
coupled. Because the solution of the coupled equations is straightforward, 
it is presented more-or-less in outline form. 

• Coupled Differential Equations (PE) 

DEf: Df r w'j" + <Mw f - w s ) - -£■ X (3) 

DE S : D sr ws" - <j> 2 (w f - w 8 ) - (4) 

where ( )' = -^-(x); tf and t s are flange and skin thicknesses; and Df r and D sr 
dx 

are flange and skin bending stiffnesses in the direction of the diagonal 
tension fold (i.e., about an inplane axis normal to the direction of the fold). 


• Boundary Condition s (BC^ at x = 0 

Note that, because the interaction between the damped waves 
emanating from the toe and heel is being neglected, boundary conditions at 
the heel are unnecessary. 


BCM): Flanae: M = 0 => w f " = 0 

(5) 

PC(2): Flange: V = 0 => 

w f ~r - (x)x=o - * Nitoe K | fr j 

2 Dfr 2 Df r \Etfr + Et S rJ 


or 


<t> ztf n 
w f = jla-L c n 

f KDfr 

(6) 

where 


r N„oeK 2 1 Etf, | 
n 2<j) z lEtfr + Etsr) 

(7) 


BCf31 : Skin: M = M t0 e = -D S rW s " 

=> W s = - (8) 

Dsr 

BCf41 : Skin: V = W toe 
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(9) 


w j"_-Wl2£. + i2kc n 

Dsr KDsr 


Solution For Wf 


Inserting 

w s = w f + J- Dfrw' f '" - x 

^Z YZ 1 

from the DE for the flange (Eq. (3)) into the DE for the skin (Eq. (4)) 
produces the following 8 th order DE for Wf: 


d 8 w f + 4p 4 d 4 w f = _ _£z 


dx' 


dx 


<!>z 2 

K 4 

Df r D sr 

L<f>z 


tfDsr + tf + ts 


Cn© 


-Kx 


where 




0z (Pfr + Dsr) 


-iL 


4DfrD 


sr 


( 10 ) 


( 11 ) 


Complementary Solution of DE (Eg, (1011 for Wf 
The solution of the homogeneous form of Eq. (10) is 
W| C) = AeP x cos (px) + BeP x sin (Px) + Ce'P x cos (Px) + De'P x sin (Px) + a + bx + cx 2 + dx 3 


as may be verified by back-substitution into the DE. 

Simplification: As described in detail above, we assume that 
the decaying waves emanating from the toe and heel die out rapidly, 
see Figure 1, so that there is no interaction between them. Thus, the 
first two terms in the above equation must vanish, and these 
conditions are satisfied by taking 


A = B = 0 


so that complementary solution for Wf simplifies to 

wj c> = Ce"P x cos (px) + De*P x sin (Px) + a + bx + cx 2 + dx 3 (12) 
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Particular Solution of DE (E g. (10)) for wf 


A particular solution of Eq. (10) is readily shown to be 


w| p) - 


K'(K , +4p 4 ) 


£ tf D sr + tf + t s 
DfrDsrL^z 


C n © 


-Kx 


Let 


m = 


1 fa 


KlK 4 +4p 4 ) D frDsr 


. tf Dsr + tf + tg 
L^z 


(dimensions of length) 


wj p) = -mC n e Kx 


(13) 

(14) 


Total Solution. Wf = Wf(°) + Wf(P) 

w f = Ce' P x cos ( Px> + De- P x sin (Px) +a + bx + cx 2 + dx 3 - mC n e Kx (15) 
• Solution For w s 


From Eq. (3): 

w s = w f + Dfr w f - x 
<l>z 0z 1 

By using Eq. (15) to eliminate Wf, and Eq. (11) to simplify the resulting 
expression, we obtain 

w s = - [Ce' P x cos ( Px) + De' P x sin ( Px] + a + bx + cx 2 + dx 3 - nC n e Kx (16) 
D$r 


where 

n = |l + K 4 m - t f 


• Determine Arbitrary Constan ts From BC’s at x = 0. 


BCm : Flange: M = 0 (Eq.(5)) => 


( 17 ) 
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(18a) 


c = p 2 D+J-mK 2 C n 


BC(2) : Flange: V = 0 (Eq. (6)) 


D = - 


Mioe + {n . m ) K 2 (Dfr+P _ sr |c n 

Dfr+Dcr) 2P 2 ' °sr i 


2p 2 (D fr+ D sr ) 


BC(3) : Skin: M = M toe (Eq. (7)) =* 

d = -lp 3 (OtO) +1 - mK 3 ]c n 

r 6 \KDfr / 


(18b) 


(18c) 


C = 


Skin: V = W, oe (Eq. (8)) 


JL 


2p 2 (Dfr+Dsr) 




+ ( _ — | j (n - m) 


'D,r+D sr / 1 2 pMP 


K - i 


iUiL 


Peel Stress At Toe (x=01 


Cn 


2p 3 K*Cfr D sr /| 

(18d) 


The peel stress at the toe is given by 
<*z = (a P eel)x=0 = <Mw f (0)-w s (0)] 


which from Eqs. (15) and (16) becomes 


g z = <|>z 


Dfr+Dsr W 
D S r / . 


+ (n - m) C n 


By substituting for C from Eq. (18d) into this equation, and performing 
some algebraic manipulation, we obtain the following expression for the 
peel stress: 


<*z = 


2p 


1 + 


D 


D fr 


[W, 


toe + 


PMtoe. + ^ztfCn 


2p 

(Dsr _ tgj 
'Dfr tfj 


1 K 4 Df r m * 

[l- K 2 + K 3 ] 

| K 

aH 

+ 

1 


<t>z . 

2p 2 2p 3 J 


(19) 


where the dimensionless parameter 
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( 20 ) 


m*=*l 
tf 

has been introduced, and m is given by Eq. (13). 

o z can be conveniently subdivided into two parts: the first, o Z f, is 
due directly to the moment Mtoe and the transverse shear Wt 0 e applied from 
the skin; and the second, <j z2 , is caused by the shear along the interface 
between the skin and flange acting at half the laminate thickness away 
from the midplane of each laminate. The two parts of the peel stress are 
given by 


<j z - G z -\ + &z2 


where 


° z1 - (TTd^dS (W|oe + 13 M|oe) 

(21) 

and 

(J Z 2 = Cn tf <>z(Fi + F 2 ) 

(22) 

and 

Fl .(l- »< 2 + K 3 |(m* Dfr K 4 . A 

\ 2 P 2 2p 3 jl ! 


P 2 P Dfr/Dsr^Dfr ■ ts/tfj 
" K \ Dfr + D sr J 
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EMANATING FROM THE HEEL DIE OUT BEFORE THEY REACH 
THE TOE (& vice versa). (ILLUSTRATED FOR PANEL Cl) 








Flange 



FIGURE 4 DIFFERENTIAL RELATIONS AND BOUNDARY CONDITIONS FOR COUPLED 
BEAM PROBLEM 





APPENDIX J 


FOUNDATION MODULI 

Expressions for the interlaminar shear stress, t xz , and peel 
stress, o z , between the skin and the attached flange are derived i n 
Appendices H and I, respectively. These expressions contain terms 
involving two “foundation moduli!” corresponding to x xz and a z . A 
simple model is used in the present appendix to obtain closed-form 
equations for these modulii. This is done by summing up the 
flexibilities of the layers between the midplanes of the skin and 
flange under the assumption that the stresses are uniform between 
the midplanes. More complex derivations are not justified within 
the framework of the overall simple approach adopted in this work. 


J.1 Foundation Modulus Corresponding to Peel Stress 


Figure 1 depicts a cross section of the attached flange and 
skin. The cross section is in the diagonal tension direction, x, and i s 
in the plane normal to the skin. A peel stress, a z , acts in the out- 
of-plane direction, z, on a unit area. Under the above assumptions, 
the extension of the midplane distance due to a z can be written as: 


5z = Gz 


ifU + 

2 'E3 /flange 


'L + _L + lf 

'E 'adhesive 'E3/face 2 'E 'core. 


In this equation, t denotes laminate thickness, E is Young’s modulus 
for the isotropic adhesive and core materials, and E3 is the out-of- 
plane modulus for the flange and skin layers. (By assuming 
transverse isotropy for tape layers, we may take E3 = E2, where E2 is 
the layer transverse modulus). 


The foundation modulus corresponding to the peel stress is 
given by 



li_L' 

2\E3)fiange 


1 


W 

M + : 


'E'adhesive 

'E3 /face < 

>Hore 


( 1 ) 


J -1 




The foundation modulus corresponding to the shear stress x X z 
is found similarly. The displacement in the x direction of the skin 
relative to the flange is 


Geff 'flange 'G'adhesive 



In this equation, G denotes the usual shear modulus for the isotropic 
adhesive and core materials, and, for the laminates 

t/Geff = X V G tti 

i th 

where Gtti is a through-the-thickness shear modulus for the i layer, 

and the summation is over the number of layers n in the laminate. In 
the case of tape layers, the values of Gtti depend on layer 
orientation. However, the value of x X z is not overly sensitive to this 
variation, and it is simplest to use G13 (subscript “1” pertains to 
the fiber direction) for all tape layers regardless of orientation. 
This simplification slightly overestimates the shear stiffness and 
so overpredicts t X z- For fabric layers G13 = G12- Thus, for both tape 
and fabric layers, we take G e ff= G13. 

The foundation modulus corresponding to the interlaminar 
shear stress is 
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J.3 Sensitivity Results 


The effect of increasing § z by 30% for the Cl panel (the most 
significantly affected panel) is to increase the maximum principal 
tensile stress at the toe, a m pt, (used to calculate failure, see 
Appendix L) by only 1.3%. Increasing <j> X z for Cl by 30% increases 
(Tmpt by 5.4% The use of the simple formulae given above is 
therefore justified. 
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SIMILAR MODEL IS USED FOR DETERMINING THE FOUNDATION MODULUS <j) X z C 
TO THE INTERLAMINAR SHEAR STRESS T xz ). 




APPENDIX K 


STRESS ANALYSIS OF THE CRITICAL PLY 


The objective of this appendix is to determine the surface tension 
stress at the toe, o x , acting in the direction of the diagonal tension fold in 
the “critical ply” of the skin. The critical ply is defined as follows: it is 

on the stiffener side of the skin; it is put into tension along the diagonal 
tension fold by the membrane load, the bending moment, and the 
transverse shear force acting in the skin at the toe of the flange; and it's 
fibers lie across the direction of the diagonal tension fold. Finally, to be 
conservative, i.e., to obtain the highest possible value for a x , we assume 
that the critical ply is at the surface immediately under the toe of the 
attached flange. The stress a x is determined from the strain components 
derived next. 

• Strains 

We define a system of inplane axes (x, y) at the middle surface of 
the skin with the x axis parallel to the direction of the diagonal tension 
fold, and the y axis being normal to it. Therefore, the strain e x in the 
direction of the diagonal tension fold, on the surface of the skin, is the 
strain normal to the fiber direction of the "critical" ply. The strain £ y 

normal to the diagonal tension fold, in the plane of the skin, on the surface 
of the skin, is the strain in the direction of the fiber of the "critical" ply. 

The internal loads in the skin laminate near the toe are a tensile 

stress resultant Nn 0 e in the x direction, a compressive stress resultant 
N2toe in the y direction, a moment Mtoe about the y axis, and a transverse 
shear force Wtoe (see Figure 1 of Appendix I). These loads cause 
membrane strains e x and e y in the midplane of the skin laminate and 

curvatures k x and K y . The surface strains in the x and y directions are 


e x = e x + -k k x 
2 

and 

e y = ey+ k Ky 


where t s is the skin thickness. 


(1) 

( 2 ) 
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From standard lamination theory [e.g., Grumman Advanced Composite 
Structures Manual], and with the assumptions that 

• the skin laminate is balanced and symmetric, 

• the bending/twisting coupling terms are small, and 

• the curvature K y about the x axis is zero, 

it is straight forward to show that curvature about the y axis is therefore 

Kx - (3) 

Dsr 

The surface strains may then be written as 

c _ Nxtoe - y N y t °Q + Mae. k_ (4) 

Et sr Et sr D sr 2 


N^ vsrl^toe. (5) 

Etsr Etsr 

In these equations, Nxtoe and Nytoe are the stress resultants in the skin 
near the toe (see Figure 1 of Appendix B). The skin laminate properties 
are relative to the diagonal tension axis system, with subscript “r” 
pertaining to the “rotated" axis system (x, y), which makes an angle of 45° 
relative to the structural axis system. Etsr is the skin membrane 
stiffness in the x direction, D S r is the neutral axis bending stiffness of 
the skin about the y axis, and v sr is Poisson’s ratio for the skin laminate. 
(v S r = v X y is defined according to v X y Ey = v yx E x> where E x and Ey are 
engineering constants for the skin laminate). 

• Stress 

The stress o x in the critical ply, in the direction of the diagonal 
tension fold, and therefore normal to the fiber direction, is related to the 
above strains by 


Ox 


Eg 

(l -V i 2 V 2 l) 


(s x + Vi 2 £y) 


( 6 ) 
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where Ei, E 2 , and v-| 2 are the in-plane Young's modulii and Poissons ratio 
for the face ply of the skin. The subscript “1" when affixed to a material 
property pertains to the fiber direction, and the Poisson’s ratio v 12 is 

defined according to v-| 2 E 2 = v 21 E-|. 

Substituting Eqs. (4) and (5) into Eq. (6), and recalling that the 
stress resultants in the skin at the toe are defined in Appendix C as 

Nxtoe = Nitoe = q (i.e., N x toe is tensile, and q>0) 


and 


Nytoe = -N2toe - -q (i.e., Nytoe is compressive, and q>0), 
we can write the equation for o x as 


o x = 


7 ^ r ((1 - V 1 2) (l +V S r) — ^ 

(1 -V 12 V 21 ) 1 Et sr 2 D sr j 


(7) 


K-3 




APPENDIX L 


FAILURE CRITERION 

Previous appendices describe a method whereby the state of stress 
in the critical ply under the toe of the stiffener flange may be rapidly, i f 
approximately, estimated. It remains to develop a failure criterion for 
this ply in order to estimate the failing strength of the stiffened panel. 

• Maximum Principal Tensile Stress 

The three stress components t xz , tfz and g x , derived in appendices H, 
I and K respectively, are combined to give the following maximum 
principal tensile stress in the surface ply immediately under the toe of 
the attached flange: 



• Failure Criterion 

Examination of photomicrographs of ply cross sections shows that 
for graphite/epoxy the placing of the fibers within such a section is 
random, and it is impossible to discern the orientation of the 
photographed section within the ply without other, external, clues. We 
may therefore assume transverse isotropy, and that the tensile strength 
at any orientation in the xz plane is equal to the transverse tensile 

strength of the layer, which can be more-or-less readily found from 
coupon tests on a laminate consisting of all 90° plies. 

Failure is assumed to occur when the maximum principal tensile 
stress, from Eq. (1), in the critical layer in the skin immediately under the 

toe reaches the transverse tensile strength F 2 U of the layer, i.e. when 

^mpt = ^2 (2) 

Failure of this layer constitutes failure of the joint, because cracking of 
the matrix allows the fibers to be pulled out of the skin surface. This 
type of failure was observed on all of the panels tested. 

The following values used in the analysis of the test panels are 
based on statistical analysis of 19 coupon test results for transverse 
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tension of IM6/3501-6 graphite/epoxy tape, in the room temperature, 
ambient, moisture ("dry”) condition [Shyprykevich, 1988]: 

[= 2 U = Mean strength = 7150 psi 
Standard Deviation = 1180 psi 
"B-basis" allowable stress = 4460 psi 

where the B-basis allowable strength is such that at least 90% of the 
transverse tensile strengths are expected to exceed the B-basis allowable 
value, with a confidence of 95%. 
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APPENDIX M 


THE SNAPPS CODE 

A Preliminary Design Code for Predicting the 

Nonlinear Response & Stiffener Separation of 
Postbuckled, Flat, Composite Shear Panels with 
Sandwich & Non-Sandwich Skins, Based on a 
Simplified Analysis Method 

The analysis methodology described previously in this report 
is based on a number of modular “strength-of-materials-type” 
models. These models are used to derive simple, closed-form 
equations that can easily be used in a "hand analysis". For 
expediency, they have been programmed into a preliminary design 
code called SNAPPS using the True Basic language. Being such a 
simple language, True Basic readily permits users to modify 
SNAPPS to suit their individual needs such as, for instance, to echo 
more input data or to provide more output. Alternatively, because 
SNAPPS requires a minimum of input data, and consists of only a 
relatively small number of executable statements (less than four 
dozen) of coding of the equations, it should be an easy task for users 
to program SNAPPS in another language or spreadsheeet (such as 
Microsoft Excel) of their choice. 

The next portion of this appendix gives a listing of the True 
Basic code. This listing also serves as a self-contained ’’User’s 
Guide” that has many annotations, along with cross-references to 
other appendices for further details. Also provided is the input data 
for three tested shear panels (A1, B1, Cl), which are the sample 
problems. These panels are described fully in Appendix A, and are 
used in Section 3 to compare predicted and measured results. 
SNAPPS output for the test panels are presented in the remaining 
portions of the current appendix. 
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M.l SNAPPS Code Listing & "User'S Guide" 


@ 000 @@ 0 @ 0000000000000000000000000000000000 @@@@@@@@@@@@@@@@@@@@@@@@@@@® 
00 SNAPPS (Version 1.0, April 1997) 00 


00 

00 

00 

00 


(Speedy Nonlnear Analysis of Postbuckled Panels in Shear) 
BY DAVE SHARP & LARRY SOBEL 

A SIMPLIFIED ANALYSIS, PRELIMINARY DESIGN CODE FOR RAPIDLY 


00 PREDICTING THE NONLINEAR RESPONSE & STIFFENER SEPARATION 


00 


OF POSTBUCKLED, FLAT, COMPOSITE SHEAR PANELS WITH SANDWICH 


00 

00 

00 

00 

00 

00 


000 000000000 00 0000000000000 0 @ 000000000000000000000000 @@@@@@@@@@@@@@@@@ 


************************************************************* 
*******************DIMENSION & DATA STATEMENTS ************** 

**********FOR GETTING RESPONSES FOR DIFF VALUES of q ******* 
************************************************************* 

The argument in the following dimension statements is the number 
of values of the shear flow, q. Response variables, such as stresses, 
are computed for each value of q. The argument can be one, 
if results for just one value of q are desired. 

DIM q(30), qratio(30), c (30) , Lb (30) , Delta(30), Mtoe(30), Wtoe(30) 

DIM Nl (30) , Lambda(30), ThetaTierodDeg (30) , DeltaThetaDeg (30) 

DIM SIGzTOE (30) , TAUxzTOE (30) , SIGxTOE(30), SIGmptTOE (30) , STRAINx (30) , STRAINy (30) 
DIM STRAINxMICRO (30) , STRAINyMICRO (30) 


************************************************************* 

******** ***********RE AD IN VALUES OF SHEAR FLOW, q ********** 
************************************************************* 

MAT READ a ! NOTE: q must always be positive (App. C) 

'DATA 634 700 791 2,800,850,900,950,1000,1050,1100,1149,1200,1300,1403,1500 
''DATA 1550,1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950,2000,2050,2100,2150,2200,2250 
'DATA 302,325,350,375,400,425,451.14,500,550,600,650,700,735,750,800,838,850 
i D ATA 900,950,1000,1050,1100,1150,1200,1250,1300,1350,1400,1450,1500 
DATA 63.5,70,75,80,85,90,95,100,150,200,250,300,350,400,467 8,500,550 

DATA 600, 650, 700, 750, 800, 850, 900, 950, 962, 1000, 1050,1100,1199 -Cl 


!A1 

!A1 

!B1 

!B1 

!C1 


******************* READ IN REST OF INPUT DATA ************** 


$####################################################################***#******* 

# COMMENT ON AXIS NOTATION (App. C) _ # 

# The long (tension) direction of the diagonal tension (DT) fold is denoted # 

# by x, and the in-plane direction normal to the DT fold is denoted by y. # 

# Most input laminate stiffness properties are referred to the DT x axis. # 

# These properties are called "rotated" properties, and are subscripted # 

§ by n r" to so designate them (rather than use double subscripts xy ) . 

# For example, Dsr is the bending stiffness of the skin IN the DT # 

# Hifpption (which is the x axis) , or eqivalently, ABOUT the y axis. * 

# ########################################## # # # ### # ### # ## ### # #################### 


READ qcr,btt, Lf 


READ tf , ts, tc, ta 

READ E3f,E3faceply,Ec,Ea 


qcr=lnitial buckling value of shear flow 
qcr is always positive (see App. C) . 
btt=Tot-to-Toe Width of Skin, measured normal 
to the direction of the stiffener 
(See Figure 1 of App. D) . 

Lf=the heel-to-toe width of the attached flange, 
measured normal to the stiffener direction. 

(See Figure 1 of App. D) . 

! Thicknesses for flange, skin, core, & adhesive. 
! Out-of-plane Young's Modulii for flange ply. 
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! face ply of skin, core, & adhesive. 

! (E3 is in dir normal to plane of skin) . 

READ G13f , G13faceply, Gc, Ga ! Out-of-plane Shear Moduli! for flange ply, 

! face ply of skin, core, & adhesive. 

! ("1" pertains to fiber direction) . 

READ Elfaceply, E2faceply, NU12faceply 

! In-plane Young's Modulii and Poisson T s Ratio for face 

! ply of skin (NU12 defined according to NU12E2=NU21El . 

! Thus, if E1>E2, then NU12 is the major Poisson ratio) . 

READ Gt , Et fr, Etsr , NUsr 

! Gt=Skin shear stiffness in the global, structural axis 

! system, x' ,y' (App. C) . This system is 

! parallel and normal to the stiffener direction. 

! Etfr=rotated membrane stiffness of flange. 

! Etsr=rotated membrane stiffness of skin. 

[ EA for a laminate is computed as follows: 

! • EA= (engineering constant for laminate)* 

! (laminate thickness) * (unit width) 

! NUsr= rotated Poisson's ratio relative to DT direction. 

! (NUsr=NUxy defined according to NUxyEy=NUyxEx) . 

READ Df r, Dsr , Dtot r 

! Dfr=rotated bending stiffness of flange 

! about neutral axis. 

! Dsr=rotated bending stiffness of skin 

! about middle surface. 

! Dtotr=rotated bending stiffness of combined (total) 

! skin and flange laminates 

! about neutral axis of total stack. 


!A1 Data, qfailtrest = 1403, qpred=1149, 
! DATA 634,9.5,0.75 
! DATA 0.072,0.1172,0.050,0.008 
! DATA 1 . 60E6, 1 . 30E6, 0 . 38E6, 0 . 30E 6 
IDATA 0.64E6,0.80E6, 0.15E6, 0.12E6 
! DATA 20.5E6, 1 .30E6, 0.35 
IDATA 0 . 2612E6, 0. 465 8E 6,0. 6302E6, 0.206 
IDATA 248.1,1215,4574 

IB1 Data, qfailtest = 838, qpred=735.0, 
IDATA 302, 9.5, 0.75 
IDATA 0.0576,0.0848,0.040,0.008 
IDATA 1 . 60E6, 1 . 30E6, 0 . 3 8E6 , 0 . 30E6 
IDATA 0 . 64E6, 0 . 80E6, 0.15E6, 0.12E6 
IDATA 20.5E6,1.30E6, 0.35 
IDATA 0.1418E6, 0.4 080E6, 0.3702E6, 0.305 
IDATA 132.5,479.9,1961 
!C1 Data, qfailtest- 962, qpred= 1199, 
DATA 63.5, 15.5, 0.75 
DATA 0.0432,0.0648,0.020,0.008 
DATA 1. 60E6, 1.30E6, 0.38E6, 0.30E6 
DATA 0.64E6,0.80E6,0.15E6,0.12E6 
DATA 20.5E6, 1.30E6,0.35 
DATA 0.1388E6, 0.3454E6, 0.3626E6, 0.306 
DATA 58.5,249, 1038 


qallow=791.2 
! qcr,btt,Lf 
tf , ts, tc, ta 
E3f ,E3f aceply, Ec,Ea 
G13f , G13faceply, Gc, Ga 
Elfaceply,E2faceply, NUl2faceply 
Gt , Etf r, Etsr, NUsr 
Df r, Dsr, Dtotr 
qallow=451 . 14 
I qcr,btt,Lf 
t f , ts, tc, ta 
E3 f , E3 f aceply , Ec, Ea 
G13f , G13faceply, Gc, Ga 
El f aceply, E2 f aceply, NUl2f aceply 
Gt , Etfr, Etsr,NUs 
Df r, Dsr, Dtotr 
qallow=4 67 . 8 
I qcr,btt,Lf 
I tf,ts,tc,ta 
! E3f,E3faceply,Ec,Ea . 

! G13f, G13faceply, Gc, Ga 
I El f aceply, E2f aceply, NU1 2 f aceply 
I Gt, Etfr, Etsr, NUsr 
! Dfr,Dsr, Dtotr 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


"&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&" 
"&&&&&&&&&&&&&&&& SNA PPS ANALYSIS OF PANEL Cl &&&&&&&&&&&&&&&" 

"&&&&&&&&&&&&&&&& (Version 1.0, April 1997) &&&£&&&&&&&&&&&" 


n************* echo selected input *************** 


"Initial Buckling Shear Flow, qcr 
"Toe-to-Toe Width of Skin, btt 
"Flange Width (normal to stiff dir.), Lf 
"Thickness of Flange, tf 


"; qcr; "lb/in. 
" ;btt; "in. " 
";Lf; "in. " 

" ;tf ; "in. " 
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PRINT "Thickness of Skin, ts 
PRINT "Thickness of Core, tc 
PRINT "Thickness of Adhesive, ta 

PRINT "Shear Stiff of Skin relative to global axes, 
PRINT "Meiab Stiff of Flange in dir of DT fold, Etfr 
PRINT "Memb Stiff of Skin in dir of DT fold, Etsr 
PRINT "Bend Stiff of Flange in dir of DT fold, Dfr 
PRINT "Bend Stiff of Skin in dir of DT fold, Dsr 
PRINT 


" ; ts ; "in. " 

" ; tc; "in . " 

" ; ta; "in . " 

Gt" ; Gt ; " lb" 

"; Etfr; "lb" 
";Etsr;"lb" 

" ;Df r; " lb-in . " 
" ; D s r ; " lb - in . " 


******** compute Geometric Variables (independent of g) ***** 


LET cosdt - 1/ sqr (2) 
LET L = btt/cosdt 

i 

LET Lfr * Lf/cosdt 


! Assumed angle of diagonal tension is45 fi 
! L - toe-to-toe length of buckle in direction of 
diagonal tension fold. 

! Lfr is width of flange in direction of 
diagonal tension fold. 


********* *get responses for diff VALUES OF q **************** 


The upper range for the index i in the following "DO" LOOP must be the same as 
the DIMENSION statement for q. The upper range can be one, 
if results for just one value of q are desired (change 
DIMENSION & MATRIX DATA statements accordingly) . 

#############################################« ######## 

# CHANGE UPPER RANGE OF INDEX FOR i # 

# TO MAKE IT EQUAL TO THE NUMBER OF q VALUES # 

####################################################### 

i 


FOR i =1 to 30 step 1 ! 


******************** STRESS RESULTANTS Nl, N2, (App. C ) ************************** 

LET qratio(i) = q(i)/qcr ! Computed for printout only 

LET Nl(i) = 2*q(i) -qcr ! Nl is the diagonal tension pull in the direction of 
1 the DT fold. 

! N2 is the compressive resultant normal to the fold. 

Recall that qcr is always positive. Thus, N2 is positive 
when it is In compression. 

******** WAVEWIDTH c & WAVELENTH Lb, (App. D) ; Max Disp, Delta (App. E) 


LET N2 - qcr 


********* 


LET c (i) = sqr ( 1 . 5) *L*sqr (N2/N1 ( i) ) / ( 1+sqr (N2/N1 (i) ) ) ! c is width of DT fold. 

LET Lb ( i) = c(i) /cosdt ! Lb is axial wavelength, computed for printout. 

LET Delta (i) = ( 2 /pi) *sqr (c (i) *L* (q (i) /Gt-qcr/Gt) ) ! Max out-of-plane disp^^^ 

I ****beNDING MOMENT Mtoe, & TRANSVERSE SHEAR FORCE AT TOE, Wtoe (App. G) 

LET Lambda (i) = sqr (Nl (i) /Dsr) 

LET ztotr = (Etfr* (ts+ta+tf/2) + Etsr*ts/2) / (Etfr+Etsr) 

! ztotr is distance of neutral axis of skin + flange model 

! from skin outer face (skin side of panel) . 

LET z = ztotr - ts/2 _ _ , . 

! z is "eccentricity" distance from middle surface of skin to 

! neutral axis of skin + flange model. , 

LET NUM = (pi*Delta ( i) / L ) * ( (Nl (i ) *Lfr A 2) / (2*Dtotr) -1 . 0) (q(i) Lfr z) /Dtotr 

LET DEN = Lambda (i) /Nl (i) + Lfr/Dtotr + pi/ (Nl (i) *L) * ( (Nl (i) *Lfr*2) / (2*Dtotr) -1 . 0) 
LET Mbar = NUM/DEN ! moment on tie rod at toe 
LET Wtoe(i)= (pi/L) * (Nl (i) *Delta (i) -Mbar) 

LET Mtoe (i) - -Mbar 


ROTATION AT TOE, Theta (for printout), (App. G) 


*********** 
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I 


LET ThetaFlangeRad - - (Mbar+q (i) * 2 ) *Lfr/Dtotr+ (Wtoe (i) *Lfr A 2) / (2*Dtotr) 

LET ThetaTierodRad - (Lambda (i) *Mbar) /N1 (i) +Wtoe (i) /N1 (i) 

LET ThetaFlangeDeg = (180/pi) *ThetaFlangeRad 
LET ThetaTierodDeg (i) = ( 180/pi) *ThetaTierodRad 

LET DeltaThetaDeg (i) - ThetaFlangeDeg - ThetaTierodDeg (i) ! check on compatibility 

l 

l ***************** INTERLAMINAR SHEAR STRESS AT TOE, TAUxzTOE (App. H) ************** 
1 

LET tface=0 . 5* (ts-tc) ! Thickness of face sheet. (This equation applies to 
! laminates without a core in which case 

! the "face" thickness is half the skin thickness) . 

LET Phixz=l/ [0.5*(tf/G13f)+ (ta/Ga) + (t face/Gl3faceply) +0.5* (tc/Gc) ] 

! Phixz is foundation modulus corresponding to shear stress (App. J) . 

LET K = sqr (Phixz* (1/Etfr+l/Etsr) ) 

LET TAUxzTOE (i) - q ( i ) *K*Et f r/ (Et f r+Et sr) ! In DT axex . 

| 

1 ********************* PEEL STRESS AT TOE, SIGzTOE (App. I) *********************** 

1 

LET Phiz=l/ [0.5* (t f/E3f ) + (ta/Ea) + (t face/E3faceply) +0.5* (tc/Ec) ] 

! Phiz is foundation modulus corresponding to Peel Stress (App. J) 

LET Beta - (Phiz* (Dsr+Df r) /4/Dsr/Dfr) A 0 . 25 

LET SIGzl * 2*Beta* (Wtoe (i) +Beta*Mtoe (i) ) / (1+Dsr/Dfr) 

! SIGzl = Peel Stress due to Mtoe & Wtoe 

LET Cn = q(i) *K A 2*Et f r/2/Phiz/ (Etsr+Etfr) 

LET mstar = (Phiz/Dsr/Dfr/K A 4 ) * (Dsr*K A 4+Phiz* (1+ts/t f ) ) / (K A 4+4*Beta A 4) 

LET FI = (l-K A 2/2/Beta A 2+K A 3/2/Beta A 3) * (mstar*K A 4*Df r/Phiz-1) 

LET F2 - 2*Beta*Dfr* (Dsr/Df r-ts/tf ) /K/ (Dfr+Dsr) 

LET SIGz2 = Cn*tf *Phiz* (F1+F2) ! Peel Stress due to shear stress. 

LET SIGzTOE (i) = SIGzl + SIGz2 ! Total Peel Stress . 

1 

! ********* STRAINS AT TOE IN CRITICAL PLY, STRAINx, STRAINx (App. K) ************* 

| 

I In what follows, (1) the x axis is in the direction of the DT fold, (2) the 
! in-plane y axis is in the direction normal to the direction of the DT fold, 

! and (3) the CRITICAL PLY is defined as follows: its fibers are normal to the 

! direction of the DT fold, it is the outermost skin ply on the stiffener side 
! of the skin, and it is conservatively assumed to be at the surface of the skin. 

LET STRAINx ( i) = q (i) /Et sr-NUsr* (-q (i) ) /Etsr+ (Mtoe (i) ) *ts/2/Dsr 

! STRAINx is the tension strain in the critical ply. 

! It is in the direction of the DT fold, 

! and is normal to the fiber direction. 

LET STRAINy ( i) * -q ( i ) /Et sr-q ( i) *NUsr /Etsr 

! STRAINy is the comp strain in the critical ply. 

! It is in the direction normal to the DT fold, 

! and it is in the fiber direction. 

LET STRAINxMICRO ( i) =STRAINx (i) *1 . 0E6 ! in units of microstrain for printout 

LET STRAINyMICRO ( i ) =STRAINy (i) *1 . 0E6 ! in units of microstrain for printout 

1 

1 *************** TENSION STRESS AT TOE IN CRITICAL PLY, SIGxTOE (App. K) *********** 

j 

LET NU21faceply = NU12faceply* (E2faceply/Elfaceply) 

LET Q22faceply = E2faceply/ (l-NU12faceply*NU21faceply) 

LET SIGxTOE (i) = Q22 faceply *STRAINx ( i ) +NU12 f aceply*Q22f aceply*STRAINy ( i) 

! SIGxTOE is the tension stress in the critical ply. 

! It is in the direction of the DT fold, 

! and is normal to the fiber direction. 

I 

I **************** m PRIN TENSION STRESS AT TOE, SIGmptTOE (App. L) *************** 

1 

LET SIGmptTOEl= (SIGzTOE ( i) +SIGxTOE ( i) ) /2 

LET SIGmptTOE 2= SQR ( ( (SIGzTOE (i) -SIGxTOE (i) ) /2) A 2+TAUxzTOE ( i) A 2) 

LET SIGmptTOE ( i ) =SIGmptTOEl+ SIGmptTOE2 

! SIGmptTOE is the max prin tension stress at the toe of the flange. 

I It acts in a plane in the DT direction that is normal to the skin. 


NEXT i 
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;**************************************★********************** 
i ******************* PRINT OUTPUT RESULTS ******************** 

I************************************************************* 

j 

! Note the following output variables (and others not shown) do not depend on q. 

! Hence, for maximum efficiency they could have been computed outside the 
! DO loop on the q values. However, this was not done because (1) we prefer to 
i compute them in the analysis sections where they belong, and (2) the code 
! runs very fast, so it really doesn't matter. 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
MAT PRINT q 

PRINT "Ratio of Applied Shear Flow to its Critical Value, q/qcr" 

PRINT " " 

MAT PRINT qratio 

PRINT "Transverse Width of Buckle, c, in." 

PRINT " 

MAT PRINT c 

PRINT "Buckle Wavelength (projected along stiffener). Lb, in." 

PRINT " 

MAT PRINT Lb 

PRINT "Maximum Out-of-Plane Displacement of Skin, Delta, in." 

PRINT " — — " 

MAT PRINT Delta 

PRINT "Axial Load in Tie-Rod model, N1 , lb/in . " 

PRINT " " 

MAT PRINT N1 

PRINT "Lambda, Wave Number in Tie-Rod Analysis, 1/in." 

PRINT " " 

MAT PRINT Lambda 

PRINT "Skin Bending Moment at Toe, Mtoe, in. -lb/in." 

PRINT " TT 

MAT PRINT Mtoe 

PRINT "Transverse Shear Force in Skin at Toe, Wtoe, lb/in." 

PRINT " 

MAT PRINT Wtoe 

PRINT "Rotation at End of Tie-rod Model, ThetaTierodDeg, deg." 

PRINT " " 

MAT PRINT ThetaTierodDeg 

PRINT "Compat . Check: Tie-rod Rotation-Skin/flange Rotation, DeltaThetaDeg, deg. 

PRINT 

MAT PRINT DeltaThetaDeg 

PRINT " In-plane Tens Strain in Critical Ply, STRAINx, micro in." 

PRINT " " 

MAT PRINT STRAINxMICRO 

PRINT " In-plane Comp Strain in Critical Ply, STRAINy, micro in." 

PRINT " 

MAT PRINT STRAINyMICRO 


W**************************************************" 

n******* SELECTED SCALAR OUTPUT RESULTS ********* 

»»***★*** (Quantities That Don't Vary with q) ******** 

»» **************************************************" 

"Toe-to-toe Length of Diagonal Tension Fold ";L;"in." 

"Thickness of Face " ; tface; " in. " 

"Thru-Thickness Foundation Modulus " ;Phiz ; "lb/sq in. /in." 

"Interlaminar Shear Foundation Modulus " ;Phixz ; "lb/sq in. /in. 1 

"Beta, Die-away Rate for Flatwise Tensile Stress" ;Beta; M/in. " 

"K, Parameter in Kuhn Shear-lag Analysis ";K;"l/in." 

«**************************************************" 

******** SELECTED MATRIX OUTPUT RESULTS ********* 

******** (Quantities That Vary with q) ********" 
„************************************************** * 


"implied Shear Flow, q, lb/in." 
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PRINT "In-plane Tens Stress in Critical Ply, SIGxTOE, Ib/sq in." 

PRINT " " 

MAT PRINT SIGxTOE 

PRINT "Peel Stress at Toe, SIGzTOE, lb/sq in." 

PRINT " " 

MAT PRINT SIGzTOE 

PRINT "Shear Lag Stress at Toe, TAUxzTOE, lb/sq in." 

PRINT " " 

MAT PRINT TAUxzTOE 

PRINT "Maximum Principal Tensile Stress Under Toe, SIGmptTOE, lb/sq in." 

PRINT " " 

MAT PRINT SIGmptTOE 


END 
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M.2 SNAPPS Output for Panel A1 


&&&&&&&&&&&&&&&& SNAPPS ANALYSIS OF PANEL A1 &&&&&&&&&&&&&&& 

&&&&&&&&&&&&&&&& (Version 1.0, April 1997) &&&&&&&&&&&&&&& 


************* ECHO SELECTED INPUT ************** 


Initial Buckling Shear Flow, qcr 

Toe-to-Toe Width of Skin, btt 

Flange Width (normal to stiff dir.), Lf 

Thickness of Flange, tf 

Thickness of Skin, ts 

Thickness of Core, tc 

Thickness of Adhesive, ta 

Shear Stiff of Skin relative to global axes, Gt 
Memb Stiff of Flange in dir of DT fold, Etfr 
Memb Stiff of Skin in dir of DT fold, Etsr 
Bend Stiff of Flange in dir of DT fold, Dfr 
Bend Stiff of Skin in dir of DT fold, Dsr 


634 lb/in. 
9.5 in. 

.75 in. 

.072 in. 
.1172 in. 

.05 in. 

.008 in. 
261200 lb 
465800 lb 
630200 lb 
248.1 lb-in. 
1215 lb-in. 


************************************************** 
******* SELECTED SCALAR OUTPUT RESULTS ******** 
******* (Quantities That Don't Vary with q) ******* 
************************************************** 


Toe-to-toe Length of Diagonal Tension Fold 
Thickness of Face 

Thru-Thickness Foundation Modulus 
Interlaminar Shear Foundation Modulus 
Beta, Die-away Rate for Flatwise Tensile Stress 
K, Parameter in Kuhn Shear-lag Analysis 


13.435 in. 

.0336 in. 

7.l0216e+6 lb/sq in. /in. 
3 . 01583e+6 lb/sq in. /in. 
9.63497 1/in. 

3.3556 1/in. 


************************************************** 
******* SELECTED MATRIX OUTPUT RESULTS ******** 
******* (Quantities That Vary with q) ******** 
************************************************** 


Applied Shear Flow, q, lb/in. 


634 

700 

900 

950 

1149 

1200 

1550 

1600 

1800 

1850 

2050 

2100 


791.2 

800 

1000 

1050 

1300 

1403 

1650 

1700 

1900 

1950 

2150 

2200 


Ratio of Applied Shear Flow to its Critical Value, q/qcr 


1 

1.41956 

1.8123 

2.44479 

2.83912 

3.23344 

Transverse Width 


1.1041 

1 . 

1.49842 

1 . 

1.89274 

2. 

2.52366 

2. 

2.91798 

2. 

3.3123 

3. 

of Buckle, 

c, in. 


24795 

1.26183 

57729 

1.65615 

05047 

2.21293 

60252 

2.68139 

99685 

3.07571 

39117 

3.47003 


850 

1100 

1500 

1750 

2000 

2250 


1.34069 

1.73502 

2.36593 

2.76025 

3.15457 

3.5489 
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8.22724 

7.83852 

7.4017 

7.36427 

7.16445 

6.98366 

6.81882 

6.66755 

6.52794 

6.39846 

6.28018 

6.16509 

5.95972 

5.77166 

5.61241 

5.53612 

5.46335 

5.39381 

5.32727 

5.2635 

5.2023 

5.1435 

5.08693 

5.03245 

4.97993 

4.92924 

4.88028 

4.83295 

4.78714 

4.74279 

Buckle Wavelength (projected along stiffener). Lb, in. 


11.6351 

11.0853 

10.4676 

10.4147 

10.1321 

9.87639 

9.64327 

9.42933 

9.23189 

9.04878 

8.88152 

8.71875 

8.42831 

8.16236 

7.93715 

7.82926 

7.72634 

7.628 

7.5339 

7.44371 

7.35717 

7.27401 

7.19401 

7.11696 

7.04269 

6.971 

6.90176 

6.83482 

6.77004 

6.70731 

Maximum Out-of-Plane Displacement of Skin, Delta 

, in. 


0 

.103849 

. 155742 

.159636 

.17961 

.196786 

.211939 

.225547 

.237929 

.249312 

.259659 

.269706 

.287649 

.304177 

.318307 

.325135 

.331689 

.337993 

.344068 

.349931 

.355599 

.361085 

.366402 

.371561 

.376574 

.381447 

.386191 

.390812 

.395318 

.399715 

Axial Load 

in Tie-Rod model, Nl, lb/in. 



634 

766 

948.4 

966 

1066 

1166 

1266 

1366 

1466 

1566 

1664 

1766 

1966 

2172 

2366 

2466 

2566 

2666 

2766 

2866 

2966 

3066 

3166 

3266 

3366 

3466 

3566 

3666 

3766 

3866 

Lambda, Wave Number in Tie 

-Rod Analysis, 1/in. 



.722365 

.794011 

.883502 

.891662 

.936678 

. 979628 

1 . 02077 

1.06032 

1.09845 

1.13529 

1.17028 

1.20561 

1.27205 

1.33703 

1.39547 

1.42465 

1.45325 

1 . 4813 

1.50882 

1.53585 

1.56242 

1.58854 

1 . 61424 

1.63953 

1.66444 

1.68899 

1.71318 

1.73703 

1.76056 

1.78379 

Skin Bending Moment at Toe 

, Mtoe, in. -lb/in. 



6.21677 

29.3103 

42.4791 

43.5241 

49.0105 

53.893 

58.3166 

62.3698 

66.112 

69.5855 

72.7592 

75.8452 

81.3274 

86.2848 

90.3948 

92.3213 

94.126 

95.8145 

97.392 

98.8631 

100.232 

101.503 

102.678 

103.763 

104.759 

105.669 

106.496 

107.242 

107.91 

108.502 

Transverse 

Shear Force in ; 

Skin at Toe, Wtoe, lb/in. 


1.4537 

25.455 

44.472 

46.237 

56.2317 

66.2565 

76.3782 

86.6284 

97.0223 

107.566 

118.048 

129.112 

151.256 

174.665 

197.243 

209.074 

221.031 

233.112 

245.313 

257.632 

270.066 

282.611 

295.266 

308.028 

320.894 

333.863 

346.932 

360.098 

373.361 

386.718 

Rotation at 

End of Tie-rod 

Model, ThetaTierodDeg, 

deg. 


-.274466 

.163239 

.419363 

.440586 

.554935 
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. 661474 

1.1328 

1.80177 

2.1918 

2.56871 


.762598 
1.22222 
1 . 88103 
2.26811 
2.64282 


.859706 

1.39315 

1.95962 

2.34393 

2.71654 


.953696 

1.56429 

2.03759 

2.41929 

2.7899 


1.04517 

1.72178 

2.11497 

2.49421 

2.86291 


Compat. Check: Tie-rod Rotation-Skin/f lange Rotation, DeltaThetaDeg, deg. 


0 

2 . 22045e-16 
-4 . 44089e-16 
8 . 88178e-16 
8 . 88178e-16 
-8 . 88178e-16 


-1 . 11022e-l6 
-6. 66134e-16 
4 . 44089e-16 
6. 66134e-16 
0 

-4 . 44089e-16 


1 . 11022e-16 
4 . 44089e-16 
4.44089e-16 
4 . 44089e-16 
0 
0 


1 . 11022e-16 
5. 55112e-16 
-2 . 22045e-16 
4 . 44089e-16 
0 

-4 . 44089e-16 


-2 . 22045e-16 
-8 . 88178e-16 
-6. 66l34e-16 
-8. 88178e-16 
4 . 44089e-16 
4 . 44089e-16 


In-plane Tens Strain in Critical Ply, STRAINx, micro in. 


1513.11 2753.22 3562 

4321.6 4630.63 4921 

5708.03 5954.46 6410 

7418.9 7601.63 7778 
8278.86 8435.82 8588 
9019.5 9155.07 9286 


89 

3630.13 

3990.42 

8 

5197.97 

5461.19 

24 

6846.45 

7230.3 

75 

7950.51 

8117.15 

22 

8736.2 

8879.92 

76 

9414.66 

9538.88 


In-plane Comp Strain in Critical Ply, STRAINy, micro in. 


-1213.27 

-1722.31 

-2198.82 

-2966.2 

-3444.62 

-3923.04 


-1339.57 

-1817.99 

-2296.41 

-3061.89 

-3540.3 

-4018 . 72 


-1514.1 

-1913.68 

-2487.78 

-3157.57 

-3635.99 

-4114.41 


-1530.94 

-2009.36 

-2684.89 

-3253.25 

-3731.67 

-4210.09 


In-plane Tens Stress in Critical Ply, SIGxTOE, lb/sq in. 


1426.08 

4872.27 

6470.24 

8359.89 

9267.2 

10018.2 


2992.93 

5233.28 

6748.36 

8555.42 

9428.98 

10151.9 


3973.71 

5570.89 

7257.75 

8743.6 

9584.77 

10280.6 


4054.08 

5888.85 

7738.87 

8924.77 

9734.77 

10404.3 


-1626.63 

-2105.05 

-2870.52 

-3348.94 

-3827.36 

-4305.78 


4482.25 

6189.83 

8156.67 

9099.21 

9879.19 

10523.2 


Peel Stress at Toe, SIGzTOE, lb/sq in. 


490.27 

1325 

2324.62 

2519 

3201.65 

3358 

4298.26 

4417 

4860.89 

4964 

5354.8 

5446 


93 

1844.36 

82 

2703.78 

27 

3648.93 

01 

4532.5 

74 

5065.97 

4 

5535.78 


Shear Lag Stress at Toe, TAUxzTOE, lb/sq in. 


904.166 998.291 1128.35 


1283.52 

1354 

1638.62 

1711 

2210.5 

2281 

2567.03 

2638 

2923.57 

2994 


82 

1426.13 

36 

1853.97 

81 

2353.11 

34 

2709.65 

87 

3066.18 


1887.05 

2878.42 

3928.59 

4644.89 

5164.67 

5623. 


1140.9 

1497.44 

2000.86 

2424.42 

2780.95 

3137.48 


2115.29 

3045.08 

4176.1 

4754.31 

5260.92 

5708.13 


1212.21 

1568.74 

2139.19 

2495.73 

2852.26 

3208.79 


Maximum Principal Tensile Stress Under Toe, SIGmptTOE, lb/sq in. 


1976.24 

5406.77 

7150.25 


3459.93 

5793.91 

7462.03 


4460.39 
6159.45 
8040 .44 


4543.99 

6506.84 

8596.52 


4992.9 

6838.56 

9088.26 
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9330.83 

9566.51 

9795.66 

10018.6 

10235.6 

10446.9 

10652.8 

10853.4 

11049 . 

11239.9 

11426 . 

11607.7 

11785 . 

11958.1 

12127.2 
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M. 3 SNAPPS Output for Panel B1 


£&£&&&&&&&&&&&&& SNAPPS ANALYSIS OF PANEL B1 &&&&&&&&&&&&&&& 

&&&&&&&&&&&&&&&& (Version 1.0, April 1997) &£&&&&&&&&&&&&& 


************* ECHO SELECTED INPUT ************** 


Initial Buckling Shear Flow, qcr 

Toe-to-Toe Width of Skin, btt 

Flange Width (normal to stiff dir.), Lf 

Thickness of Flange, tf 

Thickness of Skin, ts 

Thickness of Core, tc 

Thickness of Adhesive, ta 

Shear Stiff of Skin relative to global axes, Gt 
Memb Stiff of Flange in dir of DT fold, Etfr 
Memb Stiff of Skin in dir of DT fold, Etsr 
Bend Stiff of Flange in dir of DT fold, Dfr 
Bend Stiff of Skin in dir of DT fold, Dsr 


302 lb/in. 

9.5 in. 

.75 in. 

.0576 in. 
.0848 in. 

.04 in. 

.008 in. 
141800 lb 
408000 lb 
370200 lb 

132.5 lb-in. 
479.9 lb-in. 


************************************************** 
******* SELECTED SCALAR OUTPUT RESULTS ******** 
******* (Quantities That Don't Vary with q) ******* 
************************************************** 


Toe-to-toe Length of Diagonal Tension Fold 
Thickness of Face 

Thru-Thickness Foundation Modulus 
Interlaminar Shear Foundation Modulus 
Beta, Die-away Rate for Flatwise Tensile Stress 
K, Parameter in Kuhn Shear-lag Analysis 


13.435 in. 

.0224 in. 

8 . 73141e+6 lb/sq in. /in. 
3.663e+6 lb/sq in. /in. 
12.0413 1/in. 

4.34426 1/in. 


******* ********************************** ********* 
******* SELECTED MATRIX OUTPUT RESULTS ******** 
******* (Quantities That Vary with q) ******** 
******** ********************************** ******** 


Applied Shear Flow, q, lb/in. 


302 

325 

425 

451.14 

650 

700 

838 

850 

1050 

1100 

1300 

1350 


350 

375 

500 

550 

735 

750 

900 

950 

1150 

1200 

1400 

1450 


Ratio of Applied Shear Flow to its Critical Value, q/qcr 


1 

1.40728 

2.15232 

2.77483 

3.47682 

4.30464 

Transverse width 


1.07616 1. 
1.49384 1. 
2.31788 2. 
2.81457 2. 
3.64238 3. 
4.4702 4. 


of Buckle, c, in. 


15894 

1.24172 

65563 

1.82119 

43377 

2.48344 

98013 

3.1457 

80795 

3.97351 

63576 

4.80132 


400 

600 

800 

1000 

1250 

1500 


1.3245 

1.98675 

2.64901 

3.31126 

4.13907 

4.96689 
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8.22724 

7.93576 

7.66041 

7.41872 

7.20381 

7.01068 

6.82801 

6.52883 

6.26708 

6.0397 

5.83934 

5.66075 

5.54656 

5.50004 

5.35425 

5.252 

5.22109 

5.09874 

4.98575 

4.88092 

4.78327 

4.69197 

4.60634 

4 . 52578 

4.4498 

4.37795 

4.30986 

4.24519 

4.18367 

4.12502 

Buckle Wavelength (projected along stiffener). Lb, in. 


11.6351 

11.2229 

10.8335 

10.4917 

10.1877 

9.9146 

9.65626 

9.23315 

8.86299 

8.54143 

8.25808 

8.00551 

7.84402 

7.77823 

7.57205 

7.42745 

7.38373 

7.21071 

7.05092 

6.90266 

6.76456 

6.63545 

6.51435 

6.40042 

6.29297 

6.19136 

6.09506 

6.00361 

5.9166 

5.83367 

Maximum Out-of-Plane Displacement of Skin, Delta 

, in. 


0 

8. 37181e- 

■2 .118825 

.144207 

.164648 

.181968 

.197745 

.222798 

.244298 

.262892 

.27934 

.29413 

.303681 

.307598 

.319982 

.32878 

.331461 

. 342171 

.35222 

.361693 

.370659 

.379176 

. 387291 

.395045 

.402472 

.409602 

.41646 

.423069 

.429448 

.435615 

Axial Load 

in Tie-Rod model, Nl, lb/in. 



302 

348 

398 

448 

498 

548 

600.28 

698 

798 

898 

998 

1098 

1168 

1198 

1298 

1374 

1398 

1498 

1598 

1698 

1798 

1898 

1998 

2098 

2198 

2298 

2398 

2498 

2598 

2698 

Lambda, Wave Number in Tie 

-Rod Analysis, 1/in. 



.793283 

.851558 

. 910681 

.966192 

1.01868 

1.0686 

1.11841 

1.20601 

1.28951 

1.36793 

1.44208 

1 . 51261 

1.56008 

1.57999 

1.64461 

1.69207 

1.70678 

1.76677 

1.82479 

1.88102 

1.93562 

1.98872 

2.04043 

2.09087 

2.14012 

2.18826 

2.23537 

2.2815 

2.32672 

2.37108 

Skin Bending Moment at Toe 

, Mtoe, in.-lb/in. 



2.7567 

10.4579 

14.0668 

16.8592 

19.2241 

21.3055 

23.2567 

26.4384 

29.2109 

31.5974 

33.6624 

35.4497 

36.553 

36.9913 

38.3118 

39.1797 

39.4308 

40.364 

41.1247 

41.7243 

42.1723 

42.4772 

42.6464 

42.6866 

42.6035 

42.4026 

42.0887 

41.6661 

41.1389 

40.5107 

Transverse 

Shear Force in 

Skin at Toe, Wtoe, lb/in. 


. 644616 

9.25799 

14 . 348 

19.0492 

23.6685 

28.2997 

33.1951 

42.5468 

52.4168 

62.5919 

73.0605 

83.8079 

91.4888 

94.8189 

106.079 

114.796 

117.576 

129.296 

141.23 

153.368 

165.7 

178.219 

190.916 

203.786 

216.821 

230.017 

243.367 

256.867 

270.512 

284.298 

Rotation at 

End of Tie-rod 

Model, ThetaTierodDeg, 

deg. 


-.292593 

5.80309e-: 

2 .221353 

.352984 

.470019 
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.578464 

1.40751 

2.02249 

2.67903 

3.42151 


.685752 

1.57519 

2.06052 

2.82988 

3.56685 


.875172 

1.69058 

2.21772 

2.97947 

3.71128 


1.05896 

1.73959 

2.37309 

3.12788 

3.85486 


1.23581 

1.90124 

2.52681 

3.27521 

3.99762 


Compat. Check: Tie-rod Rotation-Skin/f lange Rotation, DeltaThetaDeg, deg. 


-1 . 11022e-l6 
1 . 11022e-16 
2 . 22045e-16 
-4 . 44089e-16 
-1 . 33227e-15 
-4 . 44089e-16 


-6. 45317e-16 
-3 . 33067e-16 
-4 . 44089e-16 
0 
0 

-8 . 88178e-16 


-1 . 38778e-16 
-8 . 88178e-16 
0 
0 

8 . 88178e-16 
8. 88178e-16 


-3 . 88578e-16 
-6 . 66134e-16 
2.22045e-16 
4 . 44089e-16 
0 
0 


-6 . 10623e-16 
6. 66134e-16 
8. 88178e-16 
4 . 44089e-16 
0 
0 


1308.15 

3380.55 

5265.46 

6415.65 

7427.37 

8329. 


Strain in 

Critical Ply, STRAINx, 

micro in. 


2069.64 

2476.62 

2811.46 

3108.53 

3645.1 

4098.44 

4519.65 

4906.76 

5599.63 

5820.49 

5912.09 

6205.02 

6480.13 

6738.84 

6982.31 

7211.53 

7630.57 

7821.77 

8001.58 

8170.5 

8477 . 52 

8616.44 

8746.12 

8866.87 

Strain in 

Critical Ply, STRAINy, 

micro in. 


-1145.66 

-1233.79 

-1321.92 

-1410.05 

-1590.32 

-1762.56 

-1938.82 

-2115.07 

-2467.59 

-2590.96 

-2643.84 

-2820.1 

-2996.35 

-3172.61 

-3348.87 

-3525.12 

-3877.63 

-4053.89 

-4230.15 

-4406.4 

-4758.91 

-4935.17 

-5111.43 

-5287.68 


-1064.59 
-1498.18 
-2291.33 
-2954.05 
-3701.38 
-4582.66 

In-plane Tens Stress in Critical Ply, SIGxTOE, lb/sq in. 


1225.72 

3742.12 

5847.97 

7051.02 

8033.86 

8811.04 


2186.24 

4046.46 

6204.96 

7116.11 

8219.26 

8924.81 


2679.04 

4561.44 

6437.76 

7374.24 

8388.95 

9025.99 


Peel Stress at Toe, SIGzTOE, lb/sq in. 


360.01 
1742.71 

2888.09 

3566.07 

4148.01 
4649. 66 


902.069 

1906.54 

3086.64 

3603.6 

4262.79 

4729.93 


1170.22 

2184.61 

3217.17 

3753.64 

4369.98 

4804.17 


Shear Lag Stress at Toe, TAUxzTOE, lb/sq in. 


687.847 

740.233 

797.174 

967.997 

1027.53 

1138.82 

1480.47 

1594.35 

1674.06 

1908.66 

1935.99 

2049.88 

2391.52 

2505.4 

2619.29 

2960.93 

3074.81 

3188.69 

Maximum Principal 

. Tensile 

Stress Under Toe 

1605.58 

2524.06 

3022.18 

4133.97 

4459.96 

5019. 

6461.36 

6875.81 

7150.29 


3077.33 

5032.48 

6533.52 

7612.4 

8543.7 

9115.06 


1385.13 

2440.4 

3271.15 

3893.97 

4469.97 
4872.61 


854.115 

1252.7 

1708.23 

2163.76 

2733.17 

3302.58 


3433.48 

5538.93 

7264.29 


3426.13 

5458.84 

6836.48 

7831.9 

8684.2 

9192.45 


1572.78 

2673.57 

3443.09 

4025.24 

4563.09 
4935.44 


911.056 

1366.58 

1822.11 

2277.64 

2847.05 

3416.46 


3798.98 

6017.35 

7629.54 
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7892.96 

7973.76 

8298.75 

8605.99 

8896.79 

9172.23 

9433.28 

9680.78 

9915.51 

10138.1 

10349.2 

10549.4 

10739.2 

10919. 

11089.2 
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M. 4 SNAPPS Output for Panel Cl 


&&&&&&&&&&&&&&«& SNAPPS ANALYSIS OF PANEL Cl &&&&&&&&&&&&&&& 

&&&&&&&&&&&&&&&& (Version 1.0, April 1997) 


************* ECHO SELECTED INPUT ************** 


Initial Buckling Shear Flow, qcr 

Toe-to-Toe width of Skin, btt 

Flange Width (normal to stiff dir.), Lf 

Thickness of Flange, tf 

Thickness of Skin, ts 

Thickness of Core, tc 

Thickness of Adhesive, ta 

Shear Stiff of Skin relative to global axes, Gt 
Memb Stiff of Flange in dir of DT fold, Etfr 
Memb Stiff of Skin in dir of DT fold, Etsr 
Bend Stiff of Flange in dir of DT fold, Dfr 
Bend Stiff of Skin in dir of DT fold, Dsr 


63.5 lb/in. 

15.5 in. 

.75 in. 
.0432 in. 
.0648 in. 
.02 in. 

.008 in. 
138800 lb 
345400 lb 
362600 lb 

58.5 lb-in. 
249 lb-in. 


************************************************** 
******* SELECTED SCALAR OUTPUT RESULTS ******** 
******* (Quantities That Don't Vary with q) ******* 
************************************************** 


Toe-to-toe Length of Diagonal Tension Fold 
Thickness of Face 

Thru-Thickness Foundation Modulus 
Interlaminar Shear Foundation Modulus 
Beta, Die-away Rate for Flatwise Tensile Stress 
K, Parameter in Kuhn Shear-lag Analysis 


21.9203 in. 

.0224 in. 

1 . 19455e+7 lb/sq in. /in. 
5 . 12601e+6 lb/sq in. /in. 
15.8456 1/in. 

5.38309 1/in. 


************************************************** 
******* SELECTED MATRIX OUTPUT RESULTS ******** 
******* (Quantities That Vary with q) ******** 
************************************************** 


Applied Shear Flow, q, lb/in. 


63.5 

70 

90 

95 

250 

300 

500 

550 

750 

800 

962 

1000 

Ratio of 

Applied Shear Flow 

i 

1.10236 

1.41732 

1.49606 

3.93701 

4.72441 

7.87402 

8.66142 

11.811 

12.5984 

15.1496 

15.748 


75 

80 

100 

150 

350 

400 

600 

650 

850 

900 

1050 

1100 

.s Critical 

Value, q/ qcr 

1.1811 

1.25984 

1.5748 

2.3622 

5.51181 

6.29921 

9.44882 

10.2362 

13.3858 

14.1732 

16.5354 

17.3228 


85 

200 

467.8 

700 

950 

1199 


I. 33858 
3.14961 
7.36693 

II. 0236 
14.9606 
18.8819 


Transverse Width of Buckle, c, in. 
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13.4234 

12.7988 

12.3881 

12.0241 

11.6977 

11.4024 

11.1332 

10.8861 

9.16312 

8.13046 

7.41248 

6.87201 

6.44424 

6.09372 

5.70489 

5.54649 

5.32658 

5.13268 

4.95991 

4.8046 

4.66393 

4.53567 

4.41808 

4.3097 

4.20938 

4.18639 

4.11614 

4.02916 

3.94777 

3.80088 

Buckle Wavelength (projected along stiffener). Lb, in. 


18.9835 

18.1003 

17.5195 

17.0046 

16.543 

16.1254 

15.7447 

15.3953 

12.9586 

11.4982 

10.4828 

9.71849 

9.11354 

8.61783 

8.06794 

7.84393 

7.53292 

7.2587 

7.01437 

6.79473 

6.59579 

6.41441 

6.2481 

6.09484 

5.95296 

5.92045 

5.8211 

5.6981 

5.58299 

5.37525 

Maximum Out 

-of-Plane Displacement of Skin, Delta, 

in . 


0 

.072971 

9 . 54 906e-2 

.112688 

.126875 

.139069 

.149821 

.159475 

.225236 

.266522 

.297461 

.322527 

.343761 

.362279 

. 384224 

.39365 

.407263 

.419822 

.431499 

.442422 

.452694 

.462398 

.471601 

.480358 

. 488717 

.490668 

.496715 

.504388 

.511764 

.525587 

Axial Load 

in Tie-Rod model 

, Nl, lb/in. 



63.5 

76.5 

86.5 

96.5 

106.5 

116.5 

126.5 

136.5 

236.5 

336.5 

436.5 

536.5 

636.5 

736.5 

872.1 

936.5 

1036.5 

1136.5 

1236.5 

1336.5 

1436.5 

1536.5 

1636.5 

1736.5 

1836.5 

1860.5 

1936.5 

2036.5 

2136.5 

2334.5 

Lambda, Wave Number in Tie- 

Rod Analysis, 1/in. 



.504995 

. 554282 

.589398 

. 622535 

. 653996 

.684011 

. 712764 

. 7404 

. 974576 

1.1625 

1.32401 

1.46786 

1.59882 

1.71983 

1.87147 

1.93934 

2.04026 

2.13641 

2.22842 

2.31678 

2.40189 

2.48409 

2.56365 

2.64081 

2.71579 

2.73348 

2.78875 

2.85985 

2.92922 

3.06194 

Skin Bending Moment at Toe, 

Mtoe, in. -lb/in. 



.288823 

1.88333 

2.4554 

2.93106 • 

3.35245 

3.73731 

4.09503 

4.43131 

7.10387 

9.06279 

10.5862 

11.7862 

12.7255 

13.4448 

14.1185 

14.3305 

14.5351 

14.5998 

14.5356 

14.3518 

14.056 

13.6552 

13.1552 

12.5611 

11.8777 

11.7008 

11.109 

10.2588 

9.33037 

7.27322 

Transverse 

Shear Force in Skin at Toe, Wtoe, Ib/in. 


4 . 13938e-2 

1.06996 

1.53571 

1.97858 

2.41703 

2.85761 

3.30313 

3.75489 

8.65248 

14.1524 

20.126 

26.4885 

33.1825 

40.1669 

50.047 

54.8888 

62.582 

70.4739 

78.5507 

86.8008 

95.214 

103.781 

112.495 

121.349 

130.335 

132.511 

139.449 

148.685 

158.04 

176.892 

Rotation at 

End of Tie-rod 

Model, ThetaTierodDeg, 

deg. 


-9 . 42543e-2 

1. 95239e-2 

5 . 86234e-2 

9 . 13684e-2 

.120798 
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.148156 

.801965 

1.65782 

2.45109 

3.09582 


.174078 
. 981229 
1.82012 
2.60509 
3.2093 


.198939 

1.15552 

1.9804 

2.75783 

3.35776 


.41893 

1.32595 

2.13888 

2.9094 

3.5053 


. 615842 

1.5521 

2.29573 

3.05986 

3.79489 


Compat. Check: Tie-rod Rotation-Skin/f lange Rotation, DeltaThetaDeg, deg^ 


0 

-3. 60822e-16 
-1 . Il022e-l6 
0 
0 

-8 . 88178e-l6 


-1 . 07553e-16 
1 . 11022e-16 
-9. 99201e-16 
-4. 44089e-16 
-4 . 44089e-16 
-4 . 44089e-l 6 


6. 93889e-18 
-4 . 44089e-16 
0 

-6.66134e-16 
-4 . 44089e-16 
8 . 88178e-16 


- 1 . 

5. 

0 

0 

4. 

0 


66533e-16 

55112e-17 


44089e-16 


1 . 52656e-16 
5 . 55112e-16 
-8 . 88178e-16 
8 . 88178e-16 
0 
0 


In-plane 

Tens 

Strain in 

Critical Ply, STRAINx, 

micro in. 


266.294 


497.183 

589.63 

669.532 

742.373 

810.459 


875.015 

936.781 

1464.62 

1899. 61 

2277 . 93 


2614.16 

2916.47 

3190.14 

3522 . 02 

3665.57 


3872.29 

4060.8 

4232.53 

4388. 69 

4530 . 3 


4658.23 

4773.26 

4876.05 

4967.21 

4987.42 


5047 . 28 

5116.73 

5176.01 

5264 . 91 

In-plane 

Comp 

Strain in 

Critical Ply, STRAINy, 

micro in. 


-228 . 712 


-252.124 

-270.132 

-288.141 

-306.15 

-324.159 


-342.168 

-360.177 

-540.265 

-720 . 353 

-900.441 


-1080.53 

-1260.62 

-1440.71 

-1684 . 91 

-1800 . 88 


-1980.97 

-2161.06 

-2341.15 

-2521 . 24 

-2701 . 32 


-2881.41 

-3061.5 

-3241.59 

-3421 . 68 

-3464.9 


-3601.77 

-3781.85 

-3961.94 

-4318 . 52 

In-plane 

Tens 

Stress in 

Critical Ply, SIGxTOE, 

lb/sq in. 



244.014 

913.198 

2571.58 

3976.74 

4696.77 

4945.53 


535.784 

989.52 

2929.52 
4164 . 99 

4781.81 

4961.2 


648.647 

1062.19 

3243.02 

4329.39 

4849.93 

4969.61 


745.076 

1671.17 

3519. 

4471.81 

4902.02 

4964.7 


832.251 

2158.5 

3841.84 

4593.83 

4938.88 

4917.66 


Peel Stress at Toe, SIGzTOE, lb/sq in. 


77.9546 

445.295 

1329.98 

2094.57 

2508.76 

2675.93 


241.615 

486.102 

1522.44 

2199.96 

2561.58 

2691.21 


303.021 
524 . 898 
1691 . 99 
2293.18 

2605.81 
2705.13 


355.08 

849.204 

1842.27 

2375.2 

2641.89 

2712.29 


401.926 

1108.97 

2019.72 

2446.83 

2670.24 

2707.55 


Shear Lag Stress at Toe, TAUxzTOE, lb/sq in. 


166.761 

183.831 

196.962 

210.092 

223.223 

236.354 

249.485 

262.616 

393.923 

525.231 

656.539 

787 . 847 

919.155 

1050.46 

1228 . 52 

1313.08 

1444 . 39 

1575.69 

1707. 

1 8 3 8 . 31 

1969. 62 

2100.92 

2232.23 

2363.54 

2494.85 

2526.36 

2626.16 

2757.46 

2888.77 

3148 . 76 

Maximum Principal Tensile 

Stress Under Toe, 

SIGmptTOE, lb/sq in. 


347 272 

624.13 

737.862 

836.719 

927.127 

1011 . 81 

1092.21 

1169.22 

1829.47 

2376 . 19 

2854.35 

3282.24 

3670.12 

4024.63 

4460.25 
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4651.14 

4929.35 

5187.27 

5426.7 

5649.13 

5855.82 

6047.87 

6226.24 

6391.76 

6545.17 

6580.26 

6687.13 

6818.24 

6939.04 

7149.65 
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